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ABSTRACT 



We study a model of weakly ionized, protostellar accretion discs that are threaded 
by a large-scale, ordered magnetic field and power a centrifugally driven wind. We 
consider the limiting case where the wind is the main repository of the excess disc an- 
Q . gular momentum and generalize the radially localized disc model of Wardle & Konigl 

?— | ' (1993), which focussed on the ambipolar diffusion regime, to other field diffusivity 

r/i \ regimes, notably Hall and Ohm. We present a general formulation of the problem for 

nearly Keplerian, vertically isothermal discs using both the conductivity-tensor and 
the multi-fluid approaches and simplify it to a normalized system of ordinary differen- 
tial equations in the vertical space coordinate. We determine the relevant parameters 
■ of the problem and investigate, using the vertical-hydrostatic-equilibrium approxima- 

V^Q ' tion and other simplifications, the parameter constraints on physically viable solutions 

Q\ , for discs in which the neutral particles are dynamically well coupled to the field already 

C*") ■ at the midplane. When the charged particles constitute a two-component ion-electron 

plasma one can identify four distinct sub-regimes in the parameter domain where the 
. Hall diffusivity dominates and three sub-regimes in the Ohm-dominated domain. Two 

iT^i of the Hall sub-regimes can be characterized as being ambipolar diffusion-like and two 

0^ . as being Ohm-like: The properties of one member of the first pair of sub-regimes are 

CO ■ identical to those of the ambipolar diffusion regime, whereas one member of the second 

pair has the same characteristics as one of the Ohm sub-regimes. All the Hall sub- 
regimes have -Brb/I-B^bl (ratio of radial-to-azimuthal magnetic field amplitudes at the 
disc surface) > 1, whereas in two Ohm sub-regimes this ratio is < 1. When the two- 
component plasma consists instead of positively and negatively charged grains of equal 
mass, the entire Hall domain and one of the Ohm sub-regimes with B r b/\B,ffo\ < 1 
disappear. All viable solutions require the midplane neutral-ion momentum exchange 
time to be shorter than the local orbital time. We also infer that vertical magnetic 
squeezing always dominates over gravitational tidal compression in this model. In a 
follow-up paper we will present exact solutions that test the results of this analysis in 
the Hall regime. 

Key words: accretion, accretion discs - MHD - ISM: jets and outflows - stars: 
formation. 



1 INTRODUCTION 

Star formation is thought to be induced by the gravitational 
collapse of the dense cores of molecular clouds (e.g. Shu et al. 
1987). During this process, angular momentum conservation 
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results in the progressive increase of the centrifugal force, a 
process that eventually halts the infalling matter and leads 
to the development of a central mass (protostar) surrounded 
by a rotationally supported, gaseous disc. In the presence of 
an angular momentum transport mechanism, mass accretion 
on to the protostar proceeds through this disc, and it is 
believed that this is how stars typically gain most of their 
mass. 
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A common feature of accreting protostellar systems is 
their association with energetic bipolar outflows that prop- 
agate along the rotation axis of the source (e.g. Bally et al. 
2007). The outflows from low-bolometric-luminosity objects 
(Lboi < 10 3 Lq) have velocities in the range ~ 150 — 400 km 
s _1 (of the order of the escape speed from the vicinity of the 
central protostar) and are typically well collimated (opening 
half-angles ~ 3 — 5 degrees on scales of 10 3 — 10 4 AU). The 
connection between accretion discs and outflows, which ex- 
tends smoothly from very low-mass stars and brown dwarfs 
(e.g. Mohanty et al. 2005) all the way up to protostars with 
masses of at least ~ 10 Mq (e.g. Corcoran & Ray 1998), is 
manifested by the observed correlations between accretion 
diagnostics (e.g. inverse P-Cygni profiles and excess emission 
in UV, IR and millimeter wavelengths) and outflow signa- 
tures (e.g. P-Cygni profiles, optical forbidden lines, thermal 
radio radiation and the presence of molecular lobes) in these 
objects (e.g. Hartigan et al. 1995). Also relevant are the cor- 
relations of the type M oc L^ ol (with S ~ 0.6 — 0.7) that 
have been inferred for both accretion and outflow rates in 
low- and high-L bo i objects (e.g. Levreault 1988) and the 
joint decline in outflow activity and disc frequency (as well 
as in the inferred accretion rate) with stellar age (e.g. Calvet 
et al. 2000, 2005). 

These features point to an underlying physical link be- 
tween accretion and outflow processes in these systems. A 
natural explanation for this connection is that outflows pro- 
vide an efficient means of transporting away the excess an- 
gular momentum and of tapping the liberated gravitational 
potential energy of the accreted matter. In this scenario, 
disc material is centrifugally accelerated by the torque ex- 
erted by a large-scale, ordered magnetic field that threads 
the disc (Blandford & Payne 1982, hereafter BP82; see also 
the reviews by Konigl & Pudritz 2000, Pudritz et al. 2007 
and Konigl & Salmeron 2009). A dynamically significant, 
ordered field component is indicated by the hourglass field 
morphology detected by polarization measurements in sev- 
eral pre-collapse molecular cloud cores on sub-parsec scales 
(e.g. Schleuning 1998; Girart et al. 2006; Kirby 2009). These 
interstellar field lines are expected to be advected inward 
when the core undergoes gravitational collapse and could 
naturally give rise to an ordered field in the resulting proto- 
stellar discs. Alternatively, an ordered field component could 
be generated locally in the disc (and star) by a dynamo pro- 
cess. The comparatively high momentum discharges inferred 
in protostellar outflow sources are consistent with this inter- 
pretation, as alternative wind driving mechanisms would fall 
far short of the indicated requirements (e.g. Konigl & Pu- 
dritz 2000). The centrifugal wind picture can also naturally 
account for the comparatively high ratio (~ 0.1) of the mass 
outflow to the mass accretion rate that is inferred in both 
the quiescent and the outburst phases of low- mass protostars 
(e.g. Ray et al. 2007; Hartmann & Kenyon 1996). 

Vertical angular-momentum transport by a large-scale, 
ordered magnetic field is expected to dominate radial 
angular-momentum transport via turbulence induced by the 
magnetorotational instability (MRI; e.g. Balbus & Hawley 
1998) when the disc is threaded by a comparatively strong 
magnetic field, corresponding to the ratio ao of the Alfven 
speed vao (where the subscript denotes the midplane) to 
the isothermal sound speed c s being not much smaller than 



unity. 1 When this condition is satisfied, angular momentum 
transport by a small-scale, disordered magnetic field is sup- 
pressed because the wavelength of the most unstable MRI 
perturbations exceeds the magnetically reduced disc scale- 
height (Wardle & Konigl e.g. 1993, hereafter WK93). On the 
other hand, when the magnetic field is comparatively weak 
(do « 1), MRI- induced radial transport should dominate. 

Previous studies have considered the aforementioned ra- 
dial and vertical angular-momentum transport mechanisms 
under a number of simplifying assumptions. The MRI has 
been analyzed in both its linear and non-linear stages (e.g. 
Sano & Stone 2002a,b; Salmeron & Wardle 2003, 2005; 
Fleming & Stone 2003; Sano et al. 2004; Desch 2004; Turner 
et al. 2007). Similarly, centrifugal wind-driving discs have 
been examined both semi-analytically, making use of self- 
similarity formulations (e.g. Wardle & Konigl 1993; Li 1995, 
1996; Ferreira 1997), and through numerical simulations 
that utilize resistive-MHD codes (e.g. Casse & Keppens 
2002; Kuwabara et al. 2005; Meliani et al. 2006; Zanni et 
al. 2007). It is, however, likely that both the radial and ver- 
tical modes of angular momentum transport play a role in 
real accretion discs. While quasi-steady disc models in which 
both mechanisms operate had been considered in the litera- 
ture (e.g. Lovelace et al. 1994; Casse & Ferreira 2000; Ogilvie 
& Livio 2001), they did not specifically link the radial trans- 
port of angular momentum with its origin in MRI-induced 
turbulence. In a first attempt in this direction, Salmeron 
et al. (2007) explored the possibility that radial and verti- 
cal transport could operate at the same radial location in 
a protostellar disc. They derived an approximate criterion 
for identifying the vertical extent of a wind-driving disc at 
the given radius that is unstable to MRI-induced turbulence 
and obtained a quantitative estimate (based on the numeri- 
cal simulations of Sano et al. 2004) of the amount of angular 
momentum that is removed radially from that region. 

An important consideration for the analysis of magnetic 
angular momentum transport in protostellar discs is that 
these discs are generally weakly ionized over most of their 
extents, so magnetic diffusivity effects are significant. This is 
a critical issue because both the centrifugal-wind and MRI- 
turbulence transport mechanisms require a minimum level 
of field-matter coupling to be effective. A complete model 
therefore needs to take into account the detailed ionization 
and conductivity structure of the disc/wind system. At low 
densities (near the disc surfaces and at large radii) the disc 
is typically in the ambipolar diffusion regime. As the den- 
sity and the shielding column increase (closer to the disc 
midplane and the central source), the fluid passes succes- 
sively through the Hall and Ohm diffusivity regimes (see 
Section 2). Note, however, that the disc parameters may 
well be such that the gas does not clearly correspond to 
any one of these limiting cases and must instead be de- 
scribed more generally using either the conductivity-tensor 
(Section 2.1) or the multi-fluid-decomposition (Section 3.14) 
formulations. Most previous treatments have simplified the 

1 Note that vertical angular-momentum transport can also 
be mediated by additional mechanisms (not considered in 
this paper), including magnetic braking (the launching of tor- 
sional Alfven waves into the ambient interstellar medium; e.g. 
Krasnopolsky & Konigl 2002), "failed" winds (which do not be- 
come supersonic or super-Alfvenic) and non-steady phenomena. 



Wind- driving protostellar accretion discs. I. Formulation and parameter constraints 3 



problem by focussing on a single conductivity regime. In par- 
ticular, WK93 and Salmeron et al. (2007) concentrated on 
the ambipolar diffusion limit (although WK93 also discussed 
Hall-current effects). WK93 were able to obtain simple an- 
alytic constraints on the parameter values for physically vi- 
able wind-driving disc solutions in this limit by modelling 
the vertical structure of the disc under the hydrostatic (negli- 
gible vertical velocity) approximation. One of their derived 
constraints yielded the minimum required degree of field- 
matter coupling, and they proceeded to verify their analysis 
by obtaining numerical solutions for discs in this strong cou- 
pling regime. 

Although it is in principle possible to seek a self-similar 
solution for the global disc/wind configuration (e.g. Konigl 
1989; Ferreira 1997), WK93 adopted a somewhat simpler 
approach in which they obtained a solution for the vertical 
structure of a radially localized (radial extent Ar < r at a 
cylindrical radius r) region of the disc and matched it to a 
self-similar (in the spherical radial coordinate R) global wind 
solution of the type originally constructed by BP82. This ap- 
proach facilitated the derivation of viability and consistency 
constraints involving exclusively the disc model parameters, 
and subsequently Li (1996) demonstrated that the basic fea- 
tures of this solution are consistent with those of a fully 
global disc/wind similarity solution. Our aim in this paper 
is to generalize the parameter constraints derived in WK93 
for the ambipolar diffusion regime to the other basic diffu- 
sivity regimes, namely Hall and Ohm. In each of the latter 
regimes we identify parameter combinations that yield dif- 
ferent sets of constraints, which correspond to solutions with 
distinct physical properties, and we use these constraints to 
delineate the Hall and Ohm sub-regimes in parameter space 
where one can derive viable, strongly coupled, wind-driving 
disc solutions. In a follow-up paper (Salmeron et al. 2009, 
hereafter Paper II) we test the results for the Hall regime 
- derived by using the hydrostatic approximation - against 
full numerical solutions obtained in the radially localized 
disc approximation of WK93. 

As noted above, the Hall and Ohm diffusivity regimes 
are likely to apply near the midplane in the inner regions 
of real protostellar discs. Typically, the surface regions from 
which the wind is driven would be sufficiently well ionized 
by cosmic rays and by the protostellar radiation field (see 
Section 2.1) to lie in the ambipolar diffusion regime (e.g. 
Salmeron & Wardle 2005; Konigl & Salmeron 2009). Nev- 
ertheless, we adopt the approximation that the entire disc 
cross section between the midplane and the base of the wind 
lies in a single diffusivity sub-regime in order to bring out 
the distinct properties of the solutions for each parameter 
range. This approximation is less appropriate for the Ohm 
regime, both because two additional regimes (ambipolar dif- 
fusion and Hall) would be encountered on the way up to 
the disc surface in this case and because the column densi- 
ties of wind-driving discs are relatively small (especially in 
comparison with those typically implied by radial-turbulent- 
transport models) and are therefore likely to contain only 
limited regions (if any at all) where the Ohm diffusivity 
dominates. We have nevertheless chosen to include the Ohm 
regime in our analysis of weakly ionized discs for complete- 
ness as well as for possible comparisons with resistive-MHD 



numerical simulations. 2 The results of this analysis could 
potentially be applicable also to the collisionally ionized in- 
nermost region of the disc, where anomalous resistivity (the 
enhanced drag between ions and electrons due to scatter- 
ing off electromagnetic waves generated by current-driven 
plasma instabilities) might develop. To obtain explicit solu- 
tions in the different diffusivity regimes, we have generalized 
the WK93 model setup by implementing a conductivity- 
tensor scheme that can be applied to any given vertical 
ionization structure. In Paper II we use a simplified ver- 
sion of this scheme to mimic the single-diffusivity conditions 
adopted in the present work. 

This paper is organized as follows. Section 2 summarizes 
the formulation, including the disc model and the governing 
equations. Section 3 describes the methodology we employ to 
reduce these equations to a set of ordinary differential equa- 
tions in the vertical coordinate z and to express them in a 
dimensionless form; it also lists the parameters of the prob- 
lem. Section 4 derives the constraints that characterize the 
parameter space where physically viable solutions exist and 
explicates them for the Hall regime, with the correspond- 
ing constraints for the Ohm regime presented in Section 5. 
The key findings of the paper are discussed in section 6 and 
summarized in Section 7. 



2 FORMULATION 
2.1 Disc model 

We model the disc as being in a steady state, geometrically 
thin, vertically isothermal, nearly Keplerian and in dynam- 
ical equilibrium in the gravitational potential of the central 
protostar. Since we expect all the physical variables to ex- 
hibit smooth spatial variations, the assumption of geometri- 
cal thinness enables us to neglect the radial derivative terms 
(|<9/<9r| ~ 1/r) in our equations in comparison with the ver- 
tical derivative terms (\d/dz\ ~ 1/h), where r and z are 
cylindrical coordinates and h(r) (<JC r) is the disc pressure 
(or density) scale-height at the radius r. 3 The disc material 
is assumed to be weakly ionized, with the abundances of 
charged species being so low that (i) the effect of ionization 
and recombination processes on the neutral gas can be ne- 
glected and (ii) the inertial, gravitational and thermal forces 
on the ionized species are negligible in comparison with the 
electromagnetic force exerted by a large-scale, ordered mag- 
netic field. Under these approximations, separate equations 
of motion for the charged species are not required, and their 
effect on the neutrals can be incorporated via a conductiv- 
ity tensor a (e.g. Cowling 1976; Wardle 1999 and references 
therein), which is a function of position {r, z}. This formula- 
tion makes it possible to systematically incorporate different 
ionized fluid components as well as the three basic field- 
matter diffusion mechanisms (ambipolar, Hall and Ohm). 



2 Note, however, that we do not continue to consider this case in 
Paper II, which concentrates on the ambipolar-diffusion and Hall 
regimes. 

3 The only exception to this approach is our retention of the 
radial derivative of the azimuthal velocity component (see 
WK93). 
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The conductivity tensor components are the Ohm, Hall and 
Pedersen terms, given by 
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respectively (e.g. Wardle & Ng 1999). In the above expres- 
sions, rij is the number density of charged species j (of mass 
rrij and total electric charge Zje), c is the speed of light, and 
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is the Hall parameter, the ratio of the gyrofrequency and the 
collision frequency of species j with the neutrals, which mea- 
sures the relative importance of the Lorentz and drag forces 
on the motion of the charged species. In these equations 

B = |B|sgn{B,} (5) 

is the signed magnetic field amplitude, with the sign intro- 
duced to keep the dependence of the Hall conductivity on 
the magnetic field polarity. 4 Also, in equation (4), 

< av >, 
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where m = 2.33 run is the mean mass of the neutral particles 
in terms of the hydrogen nucleus mass nin, p is the density 
of the neutral gas and < av >j is the rate coefficient of 
momentum exchange of species j with the neutrals. 

The above equations can be used to treat charged 
species with a wide range of masses and degrees of cou- 
pling with the neutral gas. However, to simplify the anal- 
ysis and make it more transparent, we henceforth special- 
ize to a fluid with only two charged species, one with pos- 
itive charge (which we refer to as "ions" and denote by a 
subscript i; Z\ > 0) and the other with negative charge 
("electrons," subscript e, satisfying Z a < 0). This descrip- 
tion is a good approximation both in the comparatively low- 
density regions of protostellar discs, where the two species 
can be taken as metal ions (of typical mass mi = 30 mH; 
e.g. Draine et al. 1983) and electrons, and in high-density 
(p > 1CF 13 gem -3 ) regions where the dominant charge car- 
riers are positively and negatively charged dust grains of 
equal mass (e.g. Umebayashi & Nakano 1990; Neufeld & 
Hollebach 1994; Krasnopolsky & Konigl 2002). 5 The "ion" 
and "electron" Hall parameters are related by A = q/3 e , 
where 



m e m\ + m < av > c 
mi m c + m < av >i 



(7) 



4 Note from equations (1)— (3) that, in contrast with cth, which 
depends on an odd power of B, both the Ohm and Pedersen 
conductivities arc independent of the magnetic field polarity. 

5 During the late (> 10 yr) phases of protostellar evolution, 
after the dust has largely settled to the midplane (e.g. Nakagawa 
ct al. 1981; Dullcmond & Dominik 2004; D'Alessio ct al. 2006), 
ions and electrons could dominate the charged species also in the 
dense inner regions of the disc. 



In the case of charged grains of equal and opposite charge 
as well as equal mass, q = 1, whereas in the case of ions and 
electrons we find, using 



< av > ; = 1.6 x 10 9 cm 3 s 1 (8) 



and 



<™> c ^(lxl0- 15 cm 2 )(l^y /2 (9) 

from the results of Draine et al. (1983), that q ~ 1.3 x 
10~ 4 \/r (where T denotes the gas temperature and k is 
Boltzmann's constant). In our analysis we often encounter 
the product s = f3 e /3i = f3?/q, which is sometimes referred 
to in the literature as the ion slip factor (e.g. Mitchner & 
Krugcr 1973). 

With the help of the charge-neutrality condition, 
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equations (l)-(3) become 
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respectively (e.g. Salmeron & Wardle 2003). We also make 
use of a± — \J a\ + ap , the total conductivity perpendicular 
to the magnetic field, which is given by 
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The ratios of the conductivity tensor components <th, (tp 
and ao determine the relative importance of the ambipo- 
lar, Hall and Ohm diffusion mechanisms (e.g. Wardle & Ng 
1999). These ratios, in turn, depend on the magnitudes of 
the "ion" and "electron" Hall parameters, which are func- 
tions of the neutral density and magnetic field amplitude. 
The spatial variation of the latter quantities in protostellar 
systems results in different diffusion mechanisms dominat- 
ing in different regions of the disc. We now briefly outline 
the main properties of these different regimes, although it 
must be borne in mind that some regions of the disc will 
inevitably correspond to the transition zones between two 
such regimes (or between any of the sub-regimes that we 
identify) and hence cannot be so simply classified. 

(i) Ambipolar diffusion. This regime dominates when 
ao 3> ap S> \au \ or, equivalently, when |A| 3> 1. In this case 
the magnetic field is effectively frozen into the ionized fluid 
component (the "electrons" and the "ions") and drifts with 
it through the neutrals. This regime is expected to dominate 
at relatively low densities in protostellar discs - throughout 
the vertical column in the outermost regions of the disc (e.g. 
WK93) and close to the disc surfaces at smaller radii. 

(ii) Hall diffusion (ap <g; |<th| <SC ao). These conditions 
are satisfied when the most highly mobile charged particles 
(having the highest value of \Zje/m,j\ - the electrons in the 
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case of an ion-electron plasma) are well coupled to the mag- 
netic field (\/3 c \ 3> f) even as the drift of the more massive 
charged particles of the opposite sign (i.e. the ions in an ion- 
electron plasma) through the neutrals is strongly inhibited 
by collisions (\/3i\ <§C 1). In this limit the current in an ion- 
electron plasma is dominated by the electrons, which drift 
perpendicularly to the magnetic and electric fields (thereby 
minimizing the electromagnetic force acting on them, with 
the electron-neutral drag force typically remaining negligi- 
ble). This gives rise to the Hall current, which is perpendic- 
ular to both the electric and magnetic fields. This regime 
typically dominates close to the disc midplane on "interme- 
diate" radial scales (e.g. Li 1996; Wardle 1999; Sano & Stone 
2002a; Salmeron & Wardle 2005). Note that ct h (and there- 
fore the Hall diffusivity regime) vanishes identically when 
|A| = |/?e| (i.e. wheng= 1). 

(iii) Ohm diffusion (ao w up >• |<th| or, equivalently, 
\/3 e \ <C 1). In this case all the ionized species are uncoupled 
from the field by collisions with the neutrals, so the mag- 
netic field cannot be regarded as being frozen into any of 
the fluid components. In this parameter regime one can for- 
mulate the problem in terms of a simple scalar conductivity, 
corresponding to the familiar ohmic resistivity. However, as 
we show in Section 5, one cannot in general neglect the Hall 
term in Ohm's law in classifying the physically viable disc 
solutions in this regime. The Ohm diffusivity may dominate 
in the most shielded inner regions of particularly massive 
protostellar discs (but outside the radius where collisional 
ionization of alkali metals sets in). As we noted in Section 1, 
it is also conceivable that anomalous Ohm diffusivity dom- 
inates in the collisionally ionized zone further in. In fact, 
certain models (e.g. Bell & Lin 1994; Gammie 1999) link 
the formation of a higher ionization zone in the innermost 
region of the disc to the development of FU Orionis-type 
outbursts, during which most of the accretion on to the cen- 
tral protostar may take place (e.g. Calvet et al. 2000). 

Note that the ion slip factor s = f3?/q is ^> 1 in the 
ambipolar diffusion regime and < 1 in the Ohm regime. 
This factor can be either < 1 or > 1 for a q <ti 1 two- 
component plasma in the Hall regime. This suggests that 
in the latter case one can distinguish between the ambipo- 
lar diffusion-modified Hall regime, in which the inequalities 
^fq <C \Pi \ <C 1 are satisfied, and the Ohm diffusion-modified 
Hall regime, in which the inequalities q <ti A I <S \/<7 are 
obeyed (Konigl 1997). We return to provide explicit support 
for this characterization in Section 6. 

In the fiducial minimum-mass solar nebula model 
(Hayashi 1981) the ambipolar, Hall, and Ohm diffusivity 
regimes dominate at the disc midplane on scales r > 10 AU, 
~ 1 - 10 AU and ~ 0.1 - 1AU, respectivley, for likely 
magnetic field strengths and when the effect of grains 
is neglected. If sufficiently small (< 0.1 fim) grains arc 
present at z — and the density is large enough that 
they carry a significant fraction of both the negative and 
the positive charges then the extent of the midplane Hall 
regime is strongly reduced (e.g. Wardle & Ng 1999; Wardle 
2007; Salmeron & Wardle 2008). Wind-driving accretion- 
disc models typically have lower column densities and mid- 
plane densities than either the minimum-mass model or 
standard viscous disc models (a reflection of their compara- 
tively high inward drift speeds; see Section 6): this has the 



effect of pushing both the inner and outer edges of the mid- 
plane Hall regime inward. In practice, the spatial distribu- 
tion of the different diffusivity regimes depends on the mass 
accretion rate, the dust content and the nature of the ioniza- 
tion mechanisms that affect the disc, and it could vary from 
source to source and during the evolution of any particular 
system. 



Evaluating the conductivity tensor in a real disc in- 
volves calculating the ionization balance and associated 
abundances of charged species at any given location. This 
balance is, in turn, the result of a delicate equilibrium be- 
tween ionization and recombination processes acting both in 
the gas phase and on the surfaces of dust grains (if present). 
The most relevant ionization sources outside the central 
~ 0.1 AU of the disc are nonthermal: X-rays and UV ra- 
diation, emanating from the magnetically active young star 
and its magnetosphere, as well as interstellar cosmic rays and 
radioactive decay inside the disc (e.g. Hayashi 1981; Ume- 
bayashi & Nakano 1981, 1990; Igea & Glassgold 1999). In 
the vicinity of the star the disc surface is ionized mainly 
by X-rays, but UV radiation also becomes important at 
larger distances (e.g. Semenov et al. 2004; Alexander et al. 
2005; Glassgold et al. 2005). Within the innermost ~ 0.1 AU 
the degree of ionization could be significantly raised by col- 
lisional ionization of alkali metals, particularly potassium 
(e.g. Gammie 1996; Li 1996). A detailed calculation of the 
vertical ionization and conductivity structure of a wind- 
driving disc has so far been performed only in the context 
of a radially localized treatment and without including dust 
(e.g. Konigl & Salmeron 2009) ; a generalization to a global 
disc/wind model and the incorporation of dust remain to be 
done. 



In this paper and in Paper II we adopt a simpler treat- 
ment than the one outlined above, motivated by the fact that 
our main interest is to (i) identify the regions in parameter 
space where viable wind-driving disc solutions can be found 
under weak-ionization conditions and (ii) derive the main 
properties of these solutions. As our classification scheme 
requires the specification, at the midplane of the disc, of the 
values of three distinct variables (see Section 6), we make the 
assumption that the density and ionization structure of the 
disc is such that, at the radial location under consideration, 
the parameter regime determined by the midplane values 
of the chosen three variables continues to apply throughout 
the vertical extent of the disc. In the analysis presented in 
Sections 4 and 5 of this paper we use the variables /3 e , /3 ; 
and T, the neutral-ion coupling parameter defined in equa- 
tion (82), and we assume that they remain constant with 
height at least within the quasi-hydrostatic layer. In the nu- 
merical work presented in Paper II we choose the distinct 
variables to be two independent ratios of the conductivity 
tensor components and the field-matter coupling parameter 
A (the Elsasser number; see Section 3.13), and we assume 
that they remain constant between z = and the sonic point 
above the disc surface. Notwithstanding the difference in the 
specific variables chosen in each case, these two procedures 
are effectively equivalent. 
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2.2 Governing Equations 

Our basic equations describe the conservation of mass 

dp 



dt 



+ V • (pv) = 



(15) 



and momentum 
dv 



, I — - _, JXB 

— + (v ■ V )v + - VP + V$ 

at p cp 







(16) 



for the neutral gas, as well as the evolution of the magnetic 
field 

dB 



dt 



V X E = V X (v X B) - cV X E' 



(17) 



where v is the fluid velocity and J is the current density. In 
the equation of motion (equation 16) 



PCs 



pkT 
prnn 



(18) 



is the gas pressure, c s is the isothermal sound speed and 
p = m/mH = 2.33 is the molecular weight. In addition, $ is 
the gravitational potential of the central object, given by 

GM 

(r 2 + z 2 ) 1 / 2 ' (19) 

where G is the gravitational constant and M is the mass of 
the protostar (which is treated as a point mass at the ori- 
gin of the coordinate system). Note that in equation (16) 
we have made use of the low-inertia limit for the ion- 
ized species by balancing the sum of the ion-neutral and 
electron-neutral drag forces (/ in and / cn , respectively) by 
the Lorentz force, 

J x B 



fin ~l~ f en 



(20) 



on the assumption that all other terms in the momentum 
equations for the charged particles remain small. 
In the induction equation (17) 

v X B 



E' = E + 



(21) 



is the electric field in the frame comoving with the neu- 
trals and E is the corresponding field in the inertial (labo- 
ratory) frame. The field E' is related to the current density 
J through the conductivity tensor er (see equation 24 be- 
low). Furthermore, the magnetic field satisfies the solenoidal 
condition 



V • B = 0, 



(22) 



which implies (in analogy with the inference on pv z follow- 
ing equation 35 below) that in a thin disc B z is approx- 
imately constant with height, whereas the current density 
obeys Ampere's law 



J = — V x B 

4tt 



(23) 



(where we neglected the displacement current) and Ohm's 
law, 



J = <r ■ E' ' = gqE'u + o H -B X E' ± + a P E' ± 



(24) 



Equation (24) can be inverted to yield the following 
expression for E'\ 



E' - J aH 3 X B ( — - — } 

do o\ B \ a\ g o / 



(J x B)x B 

B 2 



(25) 



The terms on the right-hand side of equation (25) are, 
from left to right, the resistive (Ohm), Hall and ambipo- 
lar diffusion contributions, respectively. Substituting equa- 
tions (11)-(14) for a two-component plasma into equa- 
tion (25) gives 

A J , A, - A J x B AA (J x B)x B 

hi 



J A - A J x B 

A + A °"c A + A cen e \Z e \ A + A Bcen c \Z c 



(26) 



where cr e is the "electron" electrical conductivity, 

a c = !0**L . (27) 
m^p 

Using A / (Bcen c | Z c |) = (c 2 7ippi) _1 in the last term of equa- 
tion (26) and taking the limit q 1 leads to 

E , = J_ + JXB { JXB)XB 

Uc cen c \Z c \ c^jipip 

which is similar in form to expressions obtained in multi- 
fluid formulations (e.g. Konigl 1989; Balbus & Terquem 
2001). 

The impact of the magnetic diffusivity on the evolution 
of the fluid can be more readily appreciated if we use the 
Ohm, Hall and ambipolar diffusivity terms, given by 

no = 7^ , (29) 



Vn = 



oh 



47T(7x OX 



and 



op 



(30) 



(31) 



47TOX Ox 

Expressing the conductivity terms in equation (25) as func- 
tions of the Ohm, Hall and ambipolar diffusivities (equa- 
tions 29-31) and substituting the resulting expression for 
E' into equation (17) leads to the following form of the in- 
duction equation (Wardle 2007): 

dB 

= V X (v x B) - V X [rjo V X B 



dt 



+ 7vh(V X B) x B + tva(V X B)x] 



(32) 



where the subscript _L again denotes the direction perpen- 
dicular to B. In the limit of a two-component plasma with 
q < 1 it can be shown (Wardle 2007) that 



and 



m = A??o 



r?A = AA»?o • 



(33) 



(34) 



Under ideal-MHD conditions all the magnetic diffusiv- 
ity terms in the induction equation tend to zero and the 
evolution of the magnetic field is fully determined by the in- 
ductive term (the first term on the right-hand side of equa- 
tions 17 and 32). 



3 METHODOLOGY 

3.1 Simplification of the equations 

We now proceed to reduce the system of equations (15) 
- (17) and (23) to a set of ordinary differential equations 
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(ODEs) in z. We closely follow the approach of WK93 and 
consider, as in that paper, a radially localized region of the 
disc that is threaded by an open magnetic field with an 
"even" symmetry (B r o = B^, = 0). Sections 3.2 - 3.8 deal 
with the mass and momentum conservation equations for 
the neutral gas and with the induction equation that de- 
scribes the evolution of the magnetic field. The remaining 
two subsections (3.9 and 3.10) are concerned with Ampere's 
and Ohm's laws, respectively. 



3.2 Continuity 

We are looking for steady-state (d/dt = 0), axisymmetric 
(d/d(j> = 0) solutions. In this limit the continuity equation 
(15) reads 



Id, . d , . 
rfr {rpVr) + d-z {pVz) 



0. 



(35) 



This expression is further simplified by neglecting the radial 
derivative term (see Section 2.1), which implies that pv z is 
constant with height. This approximation, first adopted in 
WK93, has been critiqued (e.g. Ferreira 1997) for not al- 
lowing v z to assume negative values within the disc, as it 
must do in cases (expected to be typical) in which the disc 
thickness decreases as the protostar is approached. This is- 
sue can be fully addressed only in the context of a global 
disc/wind model. However, it can be expected that any er- 
ror introduced by this approximation would be minimized if 
the upward mass flux remained small enough for v z to have 
only a weak effect on the behaviour of the other variables 
within the disc. Under the assumption that |iv| is of the 
same order of magnitude as \v$ — vk\ one can readily show 
that the condition for this to hold is that v z /c B remain <JC 1 
everywhere within the disc. This can be checked a posteri- 
ori for each derived solution of the radially localized model, 
and we do this for the solutions presented in Paper II. Note 
that the fulfillment of this condition also implies that the 
hydrostatic simplification employed in the derivation of the 
parameter constraints in Sections 4 and 5 should result in a 
good approximation to the complete numerical solution for 
the vertical structure of the disc. 



3.3 Radial component of the momentum equation 



dv r 8v r pvl dP d$ 1 

P Vr ^r + P Vz ~^ + tt + P^r = -\J$ B z - JzB^) . 

or dz r or or c 

(36) 

The dominant terms in this equation are the gravitational 
and centrifugal forces. We can thus simplify it by neglect- 
ing the dP/dr and dv r /dr terms. Because of the assumed 
geometrical thinness, we can also neglect a term of order 
(z/r) 2 v^/r in d&/dr. However, to handle the departure 
from exact Keplerian rotation, we keep the drag force (i.e. 
the Lorentz force, communicated to the neutrals by ion- 
neutral collisions; see equation 20) and the dv r /dz term. In- 
troducing these changes and using v K — « 2vk(vk — v^), 
the radial momentum equation becomes 



dv r , 2pVK , n 

pv z —, 1 (Vk - v^) 

dz r 



t>B z 



J z Ba, 



(37) 



3.4 Azimuthal component of the momentum 
equation 



dv<f, pV^Vr 

dr r 



dv<p 1 

PV Z -^- = -{J z B r 
OZ C 



(38) 



In this equation, the torque due to the ion-neutral drag (the 
right-hand side) and the inward radial transport are im- 
portant in the hydrostatic layer, whereas the vertical ma- 
terial transport becomes relevant as the wind region is ap- 
proached. It is thus appropriate to retain all the terms. Us- 
ing v$ « vk, we estimate d(rv^)/dr « «k/2. With these 
changes, equation (38) reads 



pVrVn i dv<p 1 



2r 



+ pv z —— = -(J z B r 

dz c 



JrB z 



(39) 



3.5 



Vertical component of the momentum 
equation 



dv z 



dv z 



dP <9$ 

PVr-^- + pV z — + -^-+P-^- 

dr oz oz oz 



1 



(J r B$ 



J&Br 



(40) 



The dominant terms in this equation are the thermal pres- 
sure gradient and the tidal and magnetic compression terms. 
This equation can be simplified using 



dP_ 

dz 



>.dp 

8z 



dr 



pv z 



dv z 

dz 



2 dp 
V >dz- ' ''- 



<9$ 



pv K z 



dz r r ' 

where the approximation pv z = const was used to obtain the 
first expression on the second line. The simplified equation 
is then 

(c? _^)| + ^k£ = I (J ^-J^). (41) 



3.6 Vertical component of the induction equation 



dB z 
dt 



ld_ 

r dr 



{rcE^) , 



where 



cE^ — v r B z — v z B r + cE 1 ^ 



(42) 



(43) 



In strictly steady-state discs V X E = — (l/c)dH/dt = 0, 
so E can be expressed as a gradient of a scalar function 
and hence, under the assumption of axisymmetry, £^ = 0. 
However, B z or, equivalently, the poloidal (subscript p) mag- 
netic flux function ^(r,z) (defined by 2tt^ — J J B p • dS, 
where the integral is over a spatially fixed surface that is 
threaded by the poloidal magnetic field) may still change on 
the "long" accretion time-scale r a = r/\v r \, corresponding 
to a radial drift of the poloidal magnetic field through the 
disc. The midplane radial speed of the magnetic flux surface 
at the radius of interest is given by 



VBrO = cE^jB z 



(44) 



(e.g. Nakano et al. 2002). The magnetic flux speed VBrO can- 
not be fixed arbitrarily (e.g. Ogilvie & Livio 2001) and is, 
in general, an eigenvalue of the global disc/wind problem. 
However, in the radially localized formulation of this paper 
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we write vbt-o 
Section 3.13). 



-£bc s and treat ea as a free parameter (see 



3.7 Radial component of the induction equation 



dB r 

dt 



d: 



(45) 



BE' 



dB z dv r dB r dv z 

= Vr ~^ ^ Bz ~^ v *~^ Br ~^ ^ c ~^~ 

dz dz dz oz dz 



The radial and vertical components of the induction equa- 
tion can be combined into a single relation by using the 
poloidal magnetic flux function ^ denned in Section 3.6. 
In view of the solenoidal condition on B (equation 22), vp 
satisfies 

D i a* _ i a* 

Bz — - -5— , B r — — . 

r or r oz 

Hence one can replace equations (42) and (45) by 



■at=- rcE *' 



(46) 



which implies that, if dB z /dt =fc (so that E4, =fc 0) then, 
strictly, dB r /dt need not vanish. We will, however, assume 
(as was done in WK93) that dB r /dt is identically zero. 
This approximation may be justified on the grounds that 
a nonzero value of B r can be established between the mid- 
plane (where B r = 0) and the top of the disc on the "short" 
dynamical time Td = r/v$ ~ h/c s (where h is again the disc 
scale-height) by, for example, the vertical shear in the ra- 
dial velocity field acting on B z (consider, respectively, the 
terms v z 8B r /dz and B z dv r /dz in equation 45). Although a 
fully self-consistent solution in which B r at the disc surface 
exactly matches the value determined by the global mag- 
netic field distribution outside the disc (see Section 3.13) 
may still require a contribution from the slow radial diffu- 
sion of the poloidal field, a significant change in the value 
of B r at any given radial location can potentially occur on 
a time much shorter than the accretion time. This situa- 
tion contrasts with that for the vertical field component B z 
(see equations 42 and 43) and is essentially a consequence 
of the assumed geometrical thinness of the disc. An effec- 
tively equivalent approximation was adopted in the Ogilvie 
& Livio (2001) disc model. 

Setting dBr/dt = in equation (45) implies that E^, is 
constant with height and hence (using equations 43 and 44) 
that 



CE'^ = V z B r + (vBrO - v r )B z 
at any given radius. 



(47) 



3.8 

dBj, 
dt 



Azimuthal component of the induction 
equation 



d_ 

dz 



(y (t> B z - v z Bj, 



d_ 

dr 



(VrB^ — V4.Br 



dE' r 
dz 



dE' z 
dr 



(48) 



In analogy with the case of the radial field component dis- 
cussed in Section 3.7, it can be argued that the explicit- 
time-derivative term dB^/dt in equation (48), which de- 
scribes variations in B^ over the "long" radial accretion 



time, can be neglected in comparison with the vertical advec- 
tion term v z dB^,j dz, which indicates that a measurable az- 
imuthal field component can be established over the "short" 
dynamical time, in this case primarily through the radial 
shear in the azimuthal velocity field acting on B r (the term 
Brdv^/dr in equation 48). We thus set dB^/dt — 0. In ad- 
dition, we adopt the following simplifications for the radial- 
derivative terms: 



dE z 
dr 



0, 



djVrB^) 

dr 



; 0, 



dv<t, 
dr 



' 2r 



d(v^B r 



3 B r v K 
dr 2 r 

(see Section 2.1), where we used V • B = 



dB z 
in the sec- 



ond expression. However, for consistency with the rest of 
our derivation, we continue to neglect the vertical variation 
in B z . With these approximations, and relating E' r to E r 
through equation (21), we obtain 



d_ 

dz 



(cE r 



3 B r VK 

2 r 



(49) 



3.9 Ampere's Law 

The radial and azimuthal components of Ampere's law are 
respectively, 

c dBs 



47T dz 



and 



_ c dB r 



(50) 



(51) 



where we neglected the radial derivative term (dB z /dr) in 
the last expression. The vertical component is 



J z = 



4nr 



B<t, + r 



dBj, 
dr 



• 0. 



(52) 



Our neglect of J z inside the disc is motivated by the fact 
that, under the thin-disc approximation, the magnitudes of 
the J X B terms that include this component will be small 
in comparison with the terms that involve J r and J^. 

3.10 Ohm's Law 

To simplify this equation, we begin by expressing the com- 
ponents of E' parallel and perpendicular to B as 



E'u = (E' • B)B = yB 



(53) 



and 



1 



E' ± = -^(E' xB)xB = -(yB-E'), (54) 

respectively, where y = E' ■ B/B 2 , B 2 = B 2 + B% + B 2 Z 
and B is the unit vector in the direction of B. Substituting 
equations (53) and (54) into equation (24) and using J z « 0, 
we obtain 

Jr = y(cro - <?p)B r + ^-{E'.B^ - E'^B z ) + a P E' r , (55) 
■h = y{°o - °p)B 4> + ^(E' r B z - E' z B r ) + a P E^ , (56) 
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E' = 



-B z (E' r B r + E^B^jao - op) 
Bl(a - <7 P ) + BV P 

Q-HB(E' r B c f > - E'^Br) 

B 2 z (oo-o P )+B 2 o P ' 



(57) 



3.11 System of ODEs in 2 

Summarizing, after the simplifications applied in Section 3.1, 
we have the following system of nonlinear, ordinary differ- 
ential equations in 2: 



dv r 2pV K . . 1 , 7 „ , 

/w z — 1 (vk - v$) = -(J^B Z ) . 

dz r c 



dVtf, pV r V K 1 , , „ s 



,2 2s , P»K2 1/ 7 o r R v 



J<f, 



c dB, 



<t> 



47T d2 ' 



47r dz 



3.12 Normalized Equations 



(58) 
(59) 
(60) 
(61) 
(62) 
(63) 



Equations (58)-(63) form a system of ordinary differential 
equations in v r , v^, p, E r , B$ and B r . The remaining vari- 
ables - v z , J r , J<p, E' r , E'^ and E' z - may be found alge- 
braically via equations (21) and (55)-(57), together with 
the condition pv z = const. Furthermore, as discussed above, 
under the adopted approximations B z and E^ are constant 
with height and J z w 0. 

These equations can be expressed in dimensionless form 
by normalizing the variables as follows: 



v - v K (, 



Cs 



we = 



P = 



P{r,z) 
Po(r) 



cE/B + v K r 



Cs 



cE' 

c s Bo 



(64) 



(65) 



4nh T J ... 4nh T c B a B 
3 = p , cr = , b = — , 66 

cB c z B 

where /it = c s r/vK is the tidal scale-height of the disc. Note 
that, under the adopted field symmetry, B — B z and hence 
b z — 1. The quantity w represents the normalized (by c s ) 
fluid velocity in a frame that rotates with the local Keplerian 
angular velocity 51k = vk/t. The quantity we is an analo- 
gously reduced effective flux-surface velocity. In particular, 
note that We^o = VBro/cs (see equation 44). One can simi- 
larly express the effective midplane angular velocity of the 
magnetic flux surfaces as Q.eo = —cE r o/rBo, which implies 
that WErO = —r(tt B o — f2 K )/c s . 

With these normalizations, the following dimensionless 
system of equations is obtained: 



dw r 
dz 



2 

-rj^, + 2w 4 



(67) 



dz 



1 

w z 



a . w r 
P 2 



d\np~ 

dz 



1 



1 - Wz 
dWE 



dz 



db r 
dz 

dbj, 
dz 



-j{j r b<p -j4.br) - 2 



-2 br > 



: J4> 



e' r = w Er + w,p — w z b^ , 
e'4, = -e B + w z b r - w r , 



jr = y(oo - o P )b r + -p-(e' z b^ - e'0) + a P e' r 



3<t> = y(°o - op) 



+ — (e'r 



e' z b r ) + o P e'^ , 



-{e' r b r + e'Miao ■ 



Op 1 



(ao - 5-p) + b 2 ap 



+ 



(68) 
(69) 
(70) 
(71) 

(72) 

(73) 
(74) 

(75) 
(76) 

(77) 



(cto - 5-p) + b 2 ap 

In the above expressions ao is the midplane value of a = 
va/cs, where the Alfven speed is given by 



va = 



(78) 



and £b = — WE<f, = —VBro/cs is the normalized value of the 
vertically uniform E^. 

Note in this connection that WEr, the radial component 
of the reduced electric field, is expected to be > below 
the base of the wind. This was demonstrated in WK93 for 
the ambipolar diffusion case and is also true in the Hall 
and Ohm regimes. The reason for this is that the matter 
inside the disc rotates at sub-Keplerian speeds because it 
loses angular momentum to the field, and since the motion 
of the field lines is controlled by that of the matter (given 
that the quasi-hydrostatic layer, where the bulk of the mass 
is located, is thermal pressure-dominated), f2p, must remain 
< v^/r. Hence ii)Er = —r(Q,p — S1k)/c s would be > 0. As the 
magnetic field bends away from the disc rotation axis (B r > 
above the midplane), S7k decreases with distance along 
the field line, resulting in a corresponding decline in we t - 
(Equivalently, equation 70 shows that WEr decreases with 2 
for b r > 0.) Eventually WEr vanishes at some height above 
the midplane, which roughly coincides with the location of 
the base of the wind (see Section 4.3). 

As the height of the sonic point (denoted by a subscript 
s) and the gas density at that location (or, equivalently, the 
vertical velocity at the midplane, given that pv z is constant 
with height) are not known a priori, we treat z B and w z o as 
additional variables, which satisfy 



dz s 
dz 







(79) 
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and 



dz 



0. 



(80) 



This allows us to find the position of the sonic point and the 
upward mass flux self-consistently (see Paper II) . 

3.13 Parameters 

The following parameters characterize the disc solutions in 
our model: 

(i) ao = vao/cs, the midplane ratio of the Alfven 
speed (based on the large-scale vertical field component) to 
the isothermal sound speed. This parameter measures the 
strength of the ordered magnetic field that threads the disc. 

(ii) Cs/vk = for/r, the ratio of the tidal scale-height to 
the disc radius. While this parameter, which measures the 
geometric thinness of the disc, does not appear explicitly in 
the normalized equations, it nevertheless serves to constrain 
physically viable solutions (Section 4) and is used in match- 
ing the disc solutions to the self-similar wind solutions (see 
Paper II). 

(iii) The midplane ratios of the conductivity tensor com- 
ponents: [<7p/oj_]o (or [<7h/c±]o) and [a±/ao]o- They char- 
acterize the conductivity regime of the fluid (Section 2.1). 
In general the conductivity tensor components vary with 
height, reflecting the ionization structure of the disc (e.g. 
Salmeron & Wardle 2005). However, the explicit solutions 
derived in Paper II correspond to a simplified ionization pre- 
scription wherein the above ratios remain constant with z. 

(iv) The midplane Elsasser number Ao = v^ /(QkT]±o), 
where r/± = c 2 /A-ko±_ is the "perpendicular" magnetic diffu- 
sivity. This parameter measures the degree of coupling be- 
tween the neutrals and the magnetic field, with values 3> 1 
and <C 1 corresponding to strong and weak coupling, respec- 
tively. 6 

(v) e = —Vro/c B , the normalized inward radial speed at 
the midplane. This is a free parameter of the disc solution; 
its value is determined (for given values of the other parame- 
ters) at the step where we match it to the BP82 self-similar 
global wind solution by imposing the Alfven critical-point 
constraint on the wind (see Paper II). Although the value 
of e could in principle be negative (as it in fact is in certain 
viscous disc models; e.g. Takeuchi & Lin 2002), in the con- 
text of our model formulation we generally expect e > for 
physically viable, exclusively wind-driving discs. 

(vi) £b = —vbto/cs = —cE^o/c B B z , the normalized az- 
imuthal component of the electric field E (vertically con- 
stant by equation 45) , measures the radial drift speed of the 
poloidal magnetic field lines. This parameter vanishes in a 
strictly steady-state solution but is nonzero if the magnetic 
field lines drift radially on the "long" accretion time-scale 

6 The Elsasser number A = w^/(flKf/i) is distinct from the 
Lundquist number S = vaL/t/o) and from the magnetic Reynolds 
number Rcm = VL/r/Q (where V and L are characteristic speed 
and length-scale, respectively), which have been used in similar 
contexts in the literature. This quantity was labeled by the sym- 
bol x m the MRI linear stability analyses of Wardle (1999) and 
Salmeron & Wardle (2003, 2005), where, in fact, its general form 
(which allows for variation with height) was considered. 



In a global treatment of the disc/wind problem one can 
relate eb to the value of B r at the disc surface. The latter 
is determined by the magnetic field distribution outside the 
disc, and, in the limit of a potential (V X B — 0) external 
field, can be inferred from the radial distribution of B z along 
the disc surface (e.g. Ciolek & Mouschovias 1993; Lubow et 
al. 1994; Krasnopolsky & Konigl 2002). In the numerical 
solutions presented in Paper II we simplify the calculation 
by setting eb = 0; this condition is then used to determine 
the value of B r at the disc surface. WK93, who employed a 
similar radially localized model, derived solutions for both 
positive and negative values of eb (subject to the physically 
motivated constraint eb < e) and showed that disc config- 
urations with the same value of (e — eb) were very similar 
to each other. This result suggests that setting eb equal to 
zero should not strongly impact the generality of the re- 
sults. Further discussion of this approximation is given in 
Appendix A. 



3.14 Comparison with the multi-fluid approach 

A commonly used alternative to the conductivity-tensor for- 
malism is based on writing down the separate equations of 
motion for each charged-particle species (e.g. WK93). In 
this case, instead of employing two independent ratios of 
the conductivity-tensor components, the fluid can be char- 
acterized by the electron and ion Hall parameters. Indeed, 
we already noted in Section 2 how the distinct conductiv- 
ity regimes of the fluid can be delineated directly in terms 
of p c and j3\. It is also instructive to evaluate the limiting 
forms of the magnetic coupling parameter A - the "third 
parameter" in our classification scheme of viable solutions 
(see Section 4) - in the different conductivity regimes. In 
the ambipolar diffusion limit \(3 C \ 3> \ f3\\ S> 1 and hence 



op 



ce n B 
Bft 



(81) 



(see equations 13 and 14). In this case A reduces to 

T = , (82) 



the ratio of the Keplerian rotation time to the neutral-ion 
momentum exchange time. This parameter has emerged as 
the natural measure of the field-matter coupling in the am- 
bipolar diffusion-dominated disc model of WK93 (where it 
was labelled rf), as well as in studies of the linear (e.g. Blaes 
& Balbus 1994) and the nonlinear (e.g. Mac Low et al. 1995; 
Brandenburg et al. 1995; Hawley & Stone 1998) evolution of 
the MRI in such discs. Indeed, since the ions are well cou- 
pled to the field in this limit (|/3i| > 1), the neutrals will 
be well coupled to the field if their momentum exchange 
with the charged particles (which is dominated by their in- 
teraction with the comparatively massive ions) occurs on a 
time-scale that is short in comparison with the dynamical 
time (corresponding to T > 1). 

On the other hand, in the Hall regime \/3 e \ ^> 1 and 
|/3i| <§C 1, which implies 



cr± 



(83) 



(see equations 12 and 14) and hence A = T|/3 ; |. The param- 
eter combination T|/?i|, which also figures prominently in 
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linear (e.g. Wardle 1999; Balbus & Terquem 2001) and non- 
linear (e.g. Sano & Stone 2002a,b) studies of disc MRI in the 
Hall limit, similarly has a clear physical meaning. Indeed, in 
contradistinction to the ambipolar diffusion regime, the ions 
are not well coupled to the field in this case (|/3i| < 1). In 
order for the neutrals to be well coupled to the field in the 
Hall regime it is, therefore, not sufficient for them to be well 
coupled to the ions (T > 1); rather, the product T|/3i| must 
be > 1 in this case. 

Finally, in the Ohm regime, where even the "electrons" 
are not well-coupled to the magnetic field (\/3 c \ <C 1), 



ox 



■ Op 



-0c 



(84) 



and therefore A = T/3 e /3i, implying a further tightening of 
the requirement for good coupling between the neutrals and 
the field. 

As we show in Sections 4 and 5, the parameter T plays 
a fundamental role in the theory of wind-driving discs. In 
fact, it turns out that the condition T > 1 (which, according 
to the argument given above, signifies good coupling of the 
neutrals to the field only in the ambipolar diffusion regime) 
must be satisfied by all viable configurations of this type, 
irrespective of the conductivity domain that characterizes 
the disc. We elaborate on this point in Section 6. 

One can use the definitions of the conductivity compo- 
nents to express their ratios in terms of j3 e and (3i. Special- 
izing to the case q <C 1, we have 



Op 

oh 



1 + geA 

ft 



Furthermore, in the Hall and Ohm limits we get 

ao _ ( (3 C Hall, 
oh I AT 1 Ohm. 



(85) 



(86) 



We also have 



Ao 



(87) 



which follows from the definition of the Elsasser number. 



4 PARAMETER CONSTRAINTS 

As previously shown by WK93 and Konigl (1997), viable 
wind-driving disc solutions in the strong-coupling regime 
(defined by A not being <JC 1) occupy a limited region of 
parameter space. This region is determined by the require- 
ments that (i) the flow remain sub-Keplerian within the disc, 
(ii) a wind is driven from the disc surface (i.e. a wind launch- 
ing criterion is satisfied), (iii) only the upper layers of the 
disc participate in the outflow and (iv) the rate of heating 
by Joule dissipation is bounded by the rate of gravitational 
energy release. In this section we generalize these conditions, 
which were originally derived for discs in the ambipolar dif- 
fusion regime, and apply them to discs in the Hall conduc- 
tivity domain. The corresponding constraints for discs in 
the Ohm regime are presented in Section 5. The viable so- 
lution regions can in general be delineated by the midplane 
values of two independent conductivity-tensor components 
and of the Elsasser number A (see items (iii) and (iv) in 
Section 3.13). While we employ this representation in Pa- 
per II, in this paper we use instead the parameters ip = To, 



13 = l//3io and q/3 = l//3 c o, which are appropriate to the 
two-component plasma to which we specialize and afford 
additional insights into the problem. This choice also allows 
us to make direct comparisons with the results obtained in 
WK93 for the ambipolar diffusion regime. 

The relevant constraints are derived by focussing, as in 
WK93, on the quasi-hydrostatic disc region adjacent to the 
midplane, where w\ <C 1. In this limit the vertical momen- 
tum conservation equation (69) reduces to 



dp 

dz 



alUrh - j<t>b r ) - pz. 



(88) 



The last term on the right-hand side of this equation, which 
represents tidal compression, can generally be neglected in 
comparison with the magnetic squeezing term (see Sec- 
tion 6). Applying this approximation and substituting for 
j,/, and j r from equations (71) and (72), respectively, one 
can then integrate equation (88) to deduce that 

2 

P+y(^ + ^)~ const (89) 

within the disc and, therefore, that 

b 2 rh + b% w \ (90) 

at the base of the wind (subscript b; the effective surface of 
the disc), where p C 1. Within the quasi- hydrostatic layer 
it is appropriate to set p « 1. The two other components of 
the momentum conservation relation (equations 67 and 68) 
then reduce to 



and 



2 

~ a o ■ 

w 4> « 



-2oq> ■ 



(91) 



(92) 



In the remainder of this section (as well as in Section 5) 
we assume that the ion mass density pi and the collisional 
coupling coefficient 71 (equation 6) are constant with height 
in the quasi-hydrostatic zone. Under this assumption we can 
neglect the z-variation of the quantity T within this layer 
and set it equal to To = ip. 7 



4.1 Sub-Keplerian flow below the launching 
region 

We seek to identify the regions of parameter space where 
11)4, < below the launching region, as expected for physi- 
cally viable solutions (see Section 3.12 and WK93). We begin 
by expressing j r and as functions of b r , b^, w~Er and our 
model parameters by substituting equations (21), (91) and 
(92) into equation (25). Equation (91) can then be used to 

7 In WK93 it was assumed that p; and 7; are constant from z = 
aii the way to the sonic point of the wind, and it was noted that 
this might be a reasonabic approximation to the actual conditions 
in some protostellar discs on scales > 10 2 AU, where ambipolar 
diffusion dominates. In the more general case considered in this 
paper, where we incorporate also the Hall and Ohm diffusivity 
regimes that are relevant at higher densities, a similar assumption 
would not be realistic. However, here we only make this approx- 
imation in the quasi-hydrostatic layer, where it helps to simplify 
the analysis. 
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deduce the condition that must be satisfied for the flow to 
remain sub-Keplerian within the disc. The radial and az- 
imuthal components of equation (25) are, respectively, 



crp 



jr + 

b r b<p 



b 2 ^WEr(93) 



and 



££+( 




- 




CTO / 


2 C-L 
— + 

cro 


/ crp 




\cr± 





6 2 ^b.(94) 



These relations can be cast in terms of b r , b<£, i»Er and the 
parameters if), (3 and g, which are the variables employed by 
WK93 in their analysis of the ambipolar diffusion regime. 
Specifically, inverting equations (93) and (94) and substi- 
tuting 



6 OH =&2 fM^ 
<T± \ 1 + q ) ip 

q/3 2 Ao 



fc 2 ^ = b 2 



1 + q V 



and 
we find 

/ jr 
\ U 



crp 
cr± 



cr± 



Ao 



1 + q tp 



^ (1+9) 




a 2 , if 



1"Er 



(95) 
(96) 

(97) 
(98) 



to the case q <C 1 (or, equivalently, l±g « 1), which provides 
a good approximation for discs in which grains are not the 
dominant charge carriers. It is straightforward to carry out 
this analysis in the other case of special interest, namely 
9=1, which arises when grains dominate both the "+" and 
the "-" charged components (see Section 6). In the limit 
g«l, the expression for w$ reduces to 



(105) 



-±V(2V> + (3 + b r b 4> )wv r 



(1 + q^m + b%) + (1 + /32)(1 + q*p*) + - + V) 

In the ambipolar diffusion limit (\(3 e \ 3> |/3i| 3> 1) we can 
take (3 (= 1/Ao), q/3 (= l/fto) and g/3 2 (= I/&0A0) to be 
all <JC 1. The right-hand side of equation (105) then reduces 
to the right-hand side of equation (4.3) in WK93 (where 
the minus sign in front of the numerator is a typographical 
error). On the other hand, in the Hall limit (\(3i\ < 1 < \/3 e \) 
one has /3 > 1 and q/3 <C 1, but qf3 2 may be either > 1 or 
< 1. In the Ohm regime (\/3i\ < \/3 e \ < 1) both (3 and q/3 
(and hence also q/3 2 ) are >• 1. 

Given that rasr > and K > inside the disc, the 
condition w^, < implies that 



db r 
db^ 



,2± 

jr 



2ip + f3 + b r b4, 
' l + q(3 2 + b 2 r 



(106) 



(where we used equations 98, 99 and 101 in the limit 

is < 0, or 



-brbj, <2ip + /3 = C 



(107) 



below the disc surface. The fact that b r b^ < above the 
midplane in turn implies that C must be > and hence 
that 



In these expressions 



(3 > -2-0 . 



(108) 



Al = 1 + q/3 2 + K , 
A 2 = ^(1 + q) + (3(1- q) -b r b 

A 3 = 2^(1 + q)+ (3(1 -q) + b r b 
A 4 = -(l + q/3 2 + bl) 

and 

K = (1 + qf3 2 ){b 2 r + bl) + (1 + /3 2 )(1 + q 2 (3 2 ) 



(99) 

(100) 

(101) 
(102) 

(103) 



In the remainder of this paper we assume, for simplicity, 
that ep, « (see Section 3.13 and Appendix A). Substituting 
j r from equation (92) into equation (98) and using the fact 
that WEr > (see Section 3.12) and w r < below the disc 
surface, we infer that K must be > within the disc. We 
can also substitute jj, from equation (98) into equation (91) 
to obtain an explicit expression for the sub-Keplerian flow 
condition: 

0(1 + q) [20(1 + g) + /3(1 -q)+ Kb^WEr 



2 K 
In the rest of this section and in Section 5 we also specialize 



< 0. 

(104) 



The form of equation (106) suggests that the Hall pa- 
rameter regime can be subdivided into four sub-zones, de- 
pending on how the ion slip factor so = l/q/3 2 and the El- 
sasser number Ao = ip/\(3\ compare with 1 and 1/2, respec- 
tively (Konigl 1997). In Table 1, where we list the parameter 
constraints for viable wind-driving disc solutions in the Hall 
regime, we label these sub-regions by the numerals i through 
iv. 

Combining equations (90) and (107), we infer that the 
solutions must satisfy 



,4 

Orb 



4r& 2 b + C 2 > 



(109) 



in order for the accretion flow to remain sub-Keplerian all 
the way to the top of the disc. Assuming an equality, we 
find that solutions only exist if ao > C -1 ^ 2 . This translates 
into ao > (20)~ 1//2 for Cases (i) and (ii) in Table 1 (for 
which C w 20) and a 2 , > 1/(3 for the remaining two cases 
(where C w (3). These constraints are summarized in the 
first column of inequalities in the table. 8 



8 Note that, although the ion Hall parameter can in principle 
assume a negative value, all entries of /3\q in Tables 2 and 4 that 
are not marked as absolute values are in fact > (see Section 6). 
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Table 1. Parameter constraints for wind-driving disc solutions in the limit where the Hall diffusivity dominates. Four distinct cases 
can be identified, depending on how the values of so = l/q/3 2 = /3 c oAo and of 2Ao = 2'tf>/\0\ = 2To|Aol compare with 1. The first 
inequality expresses the requirement that the disc remain sub-Kcplcrian below the wind zone (z < z b ), the second is the wind launching 
condition (the requirement that the magnetic field lines be sufficiently inclined to the vertical for centrifugal acceleration to occur), the 
third ensures that the base of the wind is located above the (magnetically reduced) density scale-height and the fourth specifies that the 
rate of Joule heating at the midplane should not exceed the rate of release of gravitational potential energy there. 



Case 


Limits 








Parameter Constraints - 


Hall Limit 








so = AioAo 


Ao 


= Tolftol 










(multi-fluid formulation) 






i 


> 1 




> 1/2 


(2T )" 1 /2 


< 


ao 
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2 


< 


eT 


< 


v K /2c s 


ii 


> 1 




< 1/2 


1/2 

PiO 


< 


ao 


< 


2(T fto) 1/2 


< 


e/20 io 


< 


T"oAo«k/cs 


iii 


< 1 




> 1/2 


(2T )-V2 


< 


ao 


< 


2 


< 


eT /3 e oAo 


< 


v K /2c s 


iv 


< 1 




< 1/2 


1/2 


< 


ao 


< 


2(T o ft ) 1/2 


< 


e/3eo/2 


< 


Tf /3io«K/c s 



Table 2. Key properties of viable disc solutions in the Hall regime. Listed, in order, are the midplane values of \db r /db^\, the magnetically 
compressed scale-height in units of the tidal scale-height (h = h/h^), the similarly normalized vertical location of the base of the wind 
z b m units of h and the normalized Joule dissipation rate j ■ e! at the midplane. 



Case 




Limits 


Solution Characteristics - Hall Limit 






so = /3 c 


Ao 


Ao = TolAol 


\db r /db<f,\ 


h 


z b /h 


U ■ e')o 


i 


> 1 




> 1/2 


2T (> 1) 


ao/eTo 


(eT ) 2 /3v / 2 


e 2 To/a 2 


ii 


> 1 




< 1/2 


Vfto (> 1) 


2a o io /t 


e 2 /6V2T o 0f Q 


e 2 /4T o 0? o a 2 


iii 


< 1 




> 1/2 


2T /3 e oAo (> 1) 


ao/eTo AoAo 


(er o e o0 io ) 2 /3V2 


t 2 T o c o0\o/a. 2 


iv 


< 1 




< 1/2 


PeO (> 1) 


2a o /e0 cO 


(e/3eo) 2 /6v / 2T Ao 


e 2 /3 c o/4T ftoag 



4.2 Wind launching criterion 

We identify the wind launching condition using a similar 
argument to the one developed in WK93. The essence of 
the argument is that, if no wind were to form, w z would 
remain small and equation (88) would be valid throughout 
the vertical extent of the disc. The right-hand side of this 
expression is < near the midplane, but if it became > 
at some height, so that 



2 

— (jrbj, ~ j<pb r 
P 



> Z, 



(110) 



this would represent an unphysical situation and would lead 
to a neccessary condition for wind launching. We obtain 
this condition by differentiating both sides of equation (110) 
with respect to z after substituting for and j r from equa- 
tion (98). In the upper layers of the disc the z derivatives of 
p, b r and 6^ can be neglected in comparison with the deriva- 
tive of WEr (given by equation 70) . The wind-launching con- 
dition then becomes (setting again p w 1) 



■ii> 2 + \W - (1 + qf) 



b 2 rh > (1 + 



? 2 )(i + g 2 /? 2 ) 



+ 7;Pip + (! + lfi 2 )b\h + Tjipq^brbb^b ■ 



(111) 



In the ambipolar diffusion regime, where (3, q/3 and q[3 2 are 
<JC 1, this expression reduces to equation (4.16) of WK93. 
That study moreover showed that viable solutions in this 
regime satisfy tp > 1 and b rb > \b$b\ (see also Section 6), 
from which it follows that this constraint further reduces to 
the ideal-MHD wind-launching criterion b rb > 1/ \/3 (BP82) 



in this case. Using b rb ~ V%/ao (see equation 90), this con- 
dition can be expressed as ao < ~ 2. The corresponding 
constraints in the Hall regime are contained in the second 
column of inequalities in Table 1. Given that the approxi- 
mation b rb w v2/ao holds in the Hall domain as well (see 
Section 4.3), it is seen that the ideal-MHD limit of the wind- 
launching criterion also characterizes Cases (i) and (iii) in 
this regime. In the other two Hall sub-regimes (Cases ii 
and iv) b rb is required to exceed (2/To/3io) 1,/2 /v / 3, which 
is > l/\/3 (i.e. the minimum inclination angle of the sur- 
face field to the rotation axis that is required for launching 
a wind is higher than the ideal-MHD value of 30°). 



4.3 Location of the base of the wind 

In light of both theoretical and observational arguments (e.g. 
Konigl & Pudritz 2000), only a small fraction of the disc 
material should participate in the outflow. This condition is 
implemented by requiring that the base of the wind (which 
we identify with the height above which the azimuthal ve- 
locity becomes super-Keplerian) be located above the mag- 
netically reduced density scale-height (i.e. z b /h > 1). When 
this condition is violated, B^, changes sign within the disc, 
reflecting the attempt by the field to transfer angular mo- 
mentum to matter before the outflow is initiated (see Fig. 6 
in WK93 and Paper II). 

To formulate this constraint, we begin by estimating 
the vertical location of the base of the wind. In the hydro- 
static approximation, the azimuthal velocity increases to a 
Keplerian value (i.e. — 0) when ii)E r decreases to (see 
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equations 67 and 98). By equation (70), this occurs at 

2 IDErO 
2b « 77- 



3 & r b 



(112) 



Next, we substitute j r0 from equation (92) and w r0 = — e 
into equation (98) to get 

e ^ 2 + (5/2)^/3 + (1 + /3 2 )(1 + g 2 /3 2 ) 

WEr0 = W 1TW 2 ' ( } 

where, in the Hall regime (/3 2 ^> 1 and q 2 p 2 <JC 1), the 
product (1 + /3 2 )(1 + q 2 P 2 ) on the right-hand side reduces to 
/3 2 . 9 As an aside, we note that equation (113) can be used 
to derive a parameter consistency constraint based on the 
requirement that w>ErO (or, equivalently, the midplane value 
of K; see equation 103) must be > (see Section 3.12). 
Indeed, treating this expression as a quadratic equation for 
P and setting it equal to 0, we find that viable solutions in 
the Hall regime require (3 < —2tp or P > —tp/2. In view of 
equation (108), it follows that 



P > —tp/2 in Hall regime. 



(114) 



In the wind-driving discs that are of interest to us the 
density scale-height (h) is generally reduced by magnetic 
squeezing (h = h/fiT < 1; see Section 6 and WK93). To 
obtain an expression for h we first approximate the sur- 
face values of b r and 6^ by b r b ~ Zb(db r /dz) and b^b « 
Zb(db^/dz)o, respectively. 10 Then, using equation (106), we 
deduce 



(115) 



brb 




db r 


2tP + /3 


bif,b 




db^ 


o 1 + qP 2 



which is readily verified to be > 1 in all the Hall sub-regimes. 
These results are presented in column 4 of Table 2, which 
summarizes the key properties of our wind solutions for the 
various Hall sub-regimes considered in this section. Neglect- 
ing the b^ term in equation (88) and approximating p w 1 
and b r « z(db r /dz)o, we find 



dp 

dz 



db r 

dz 



+ 1 



(116) 



As we show in Section 6, magnetic compression (represented 
by the first term on the right-hand side of equation 116) 
dominates tidal squeezing (the second term on the right- 
hand side) in all viable wind-driving disc solutions. Defining 
the scale-height h by relating it to an effective Gaussian 
density distribution, so that 

dp z 
d2 ~ ~& ' 

we thus infer 



(117) 



fflo 



db r 
dz 



(118) 



By substituting for (db r /dz)o in equation (118) from 
equations (71) and (98), we can express h in terms of the 
disc parameters as 



h 



2a l + qf3 2 



2t(j+P 



(119) 



Similarly, combining equations (90), (112), (113) and (119) 
we get the desired expression for the location of the base of 
the wind in the Hall regime: 



Zb 

h 



e{2tp + P) tp 2 + (5/2)V>/3 + f3 2 



&V2tp 



(1 + qP 2 ) 2 



(120) 



Expressions for h valid for Cases (i) — (iv) are listed 
in column 5 of Table 2. Equation (120), which is listed in 
column 6 of that table, can be used to express the con- 
straint Zb/h > 1 in each of the four Hall sub-regimes. This 
constraint is shown in the third column of inequalities in 
Table 1. 



4.4 Dissipation rate 

This constraint, represented by the last set of inequalities 
in Table 1, limits the rate of heating by Joule dissipation at 
the midplane to less than the rate of gravitational potential 
energy released at that location (Konigl 1997), 

e Vk 
2a 2 c s ' 

where, by equations (26), (27) and (52), 



U ■ e')o < 



(121) 



(i-e')o = #(l + 4/3 2 )OV 2 o+^o) 



tp 



^Pal 



(1 + qP 2 



1 + 



2tP + P 
1 + qP 2 



(122) 



The expression (122) is shown in column 7 of Table 2 for 
each of the Hall sub-regimes. 



5 THE OHM REGIME 

In this section we consider the parameter constraints on vi- 
able wind-driving disc solutions for discs in the Ohm diffu- 
sivity regime. The corresponding constraints for discs in the 
Hall regime were derived in Section 4 using the hydrostatic 
approximation, and here we either directly apply the re- 
sults presented in that section or else further generalize them 
to cover also the Ohm domain. We again adopt the multi- 
fluid formalism in a weakly ionized "ion"- "electron" gas with 
q <§C 1, in which the Ohm regime is defined by the inequali- 
ties \p i0 \ < \Pco \ < 1 (or, equivalently, \0\ > q\P\ > 1). We 
also continue to assume that ee ~ 0. 



9 A more general expression for uiErO; m which no restriction is 
placed on the value of q, is given by equation (All) in WK93. 

10 We will not use these expressions to estimate the actual mag- 
nitudes of the transverse magnetic field components at the base 
of the wind (rather than just their ratio) both because the as- 
sumption of a linear scaling of the field components with z is not 
expected to remain accurate above a density scale height and be- 
cause of the approximation involved in the definition (112) of Zb\ 
instead, we will use equation (90) for this purpose. 



5.1 Sub-Keplerian flow below the launching 
region 

All the results derived in Section 4.1 are applicable in this 
regime. In particular, a necessary condition for a solution to 
exist is that the inequality ao > C -1 ^ 2 is satisfied, where C, 
given by equation (107), is equal to 2tp (where tp = To) or to 
P in the limits where tp/\P\ is 3> 1/2 or <g 1/2, respectively. 
The magnitude of this parameter combination thus provides 
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a natural classification criterion, similarly to the situation in 
the Hall regime. However, in contrast with the Hall case, in 
which the factor q(3 2 in the denominator of equation (106) 
can be either > 1 or < 1, in the Ohm limit q/3 ^> 1 and 
this parameter combination cannot be used to identify the 
relevant parameter sub-regimes. A suitable combination can 
nevertheless be found by using the relation (106) to deter- 
mine which transverse magnetic field component (radial or 
azimuthal) dominates at the disc surface (see equation 115). 
In the limit \f3\ S> 2t[> we have l&rb/fr^bl = l/<7/3 <C 1, 
but in the opposite limit (2ip 3> |/3|) this ratio is given by 
\b r b/b(f,b\ — 2ip/q/3 2 , which can be either > 1 or < 1. This 
suggests that the parameter combination tp/q/3 2 could serve 
to sub-divide the parameter regime i/>//3 > 1/2. As we verify 
in the subsequent subsections, this is indeed a proper choice. 
This result is not surprising given the fact that, in the Ohm 
regime, 4>/qf3 2 is equal to the Elsasser number Ao (defined 
in Section 3.13) and that Ao plays a similar role in the Hall 
regime (where it is equal to tp/\/3\). 



5.2 Wind launching criterion 

In analyzing the inequality given by equation (111) in the 
Ohm limit, we consider separately the cases "i/>/|/3| > 1/2 and 
ip/\(3\ < 1/2. As we noted in Section 5.1, the first case can 
be sub-divided according to whether Ao is > 1/2 or < 1/2. 
When Ao > 1/2, the dominant terms on both the left-hand 
side and the right-hand side of equation (111) are the ones 
containing the factor ip 2 , and we recover the BP82 launch- 
ing condition of a centrifugally driven, ideal-MHD wind from 
a Keplerian disc, 6 rb > I/a/3- However, the application of 
the above equation to the second sub-regime, in which the 
inequalities < 2ip < q/3 2 are satisfied, is less straight- 
forward. Indeed, while one can readily verify that the term 
involving ip again dominates the left-hand side of equa- 
tion (111) in this case, on the right-hand side of this equa- 
tion either the first term or the last term could potentially 
dominate. Assuming that the first term (~ q 2 P i ) in fact 
dominates, one can substitute b 2 h ~ {2/a 2 ] )(2tp / 'q[3 2 ) 2 from 
equations (90) and (115) into equation (111) to infer that 
the factor Qip 2 / 'a\q 2 j3 A is > (qf3 2 /2tp) 2 > 1. But in the cho- 
sen parameter regime this factor is equal to the ratio of the 
absolute magnitude of the last term on the right-hand side 
of equation (111) to the first term on the right-hand side 
of this equation, which is a contradiction (as we assumed 
that the first term on the right-hand side dominates). A 
self-consistent solution can thus be obtained only if the ra- 
tio of the absolute values of the last and first terms on the 
right-hand side of equation (111) is required from the start 
to exceed 1, in which case the inequality expressed by this 
equation would be trivially satisfied (with the left-hand side 
being > and the right-hand side - which is proportional 
to b rb b<j, b - being < 0). This requirement, in turn, provides 
an effective wind- launching condition on the inclination of 
surface poloidal field: b r h > 2/\/3. 11 

Turning now to the case ip/\/3\ < 1/2, we find that it 
is again possible to choose between two alternatives: either 



ip > q\/3\ or ip < q\/3\. The first alternative is similar to 
the just-discussed ip/\f3\ > 1/2, Ao < 1/2 sub-regime in that 
one can immediately identify the dominant term on the left- 
hand side of equation (111) (which in this instance involves 
the factor Zip 13/ 2) even as the choice between the first and 
last terms on the right-hand side of this equation appears to 
be unclear. The assumption that the first term on the right- 
hand side dominates again leads to a contradiction, so in this 
case, too, one is led to impose the requirement that the last 
term always dominates the first term on the right-hand side 
of equation (111). This requirement, in turn, again implies 
an effective wind-launching criterion, which in this case is 
brb > [(4/3)(/3/2t/>)] 1/2 (> 2/V3). 12 The second alternative, 
corresponding to the parameter regime q\(5\ < ip < \f}\/2, 
does not allow the launching condition to be satisfied since 
it implies that, while the right-hand side of equation (111) 
remains > 0, the left-hand side is dominated by the last 
term and is therefore < 0. We therefore exclude this pa- 
rameter regime from further consideration in the ensuing 
subsections, in which we continue to verify that the three 
other Ohm sub-regimes that we have identified are com- 
patible with wind-driving disc solutions that obey all the 
relevant physical requirements. 



5.3 Location of the base of the wind 

In the Hall regime we deduced that all viable solutions sat- 
isfy l&rb/fr^bl ~ \db r j 'd&^lo 3> 1, and we used this approxima- 
tion in deriving the expressions (118) and (120) for the nor- 
malized density scale height h and disc height Zb. However, 
as already noted in Section 5.1, in the Ohm case this inequal- 
ity could in principle be reversed, in which case the above 
expressions would need to be modified. This can be done 
through a straightforward generalization of the derivations 
presented in Section 4.3. Using equation (115) and defining 



D(1>,/3) 



fdbr\ 2 = (2 



7 2</> + P x 2 



+ q/3 2 



we obtain 



h= 2 -^[l + D^,P)]-^, 
brb ^^.[ 1 + 1 /D(xb,/3)}- 1/2 



(123) 



(124) 



(125) 



and 



5 b „ <? V 2 + (5/2)^/3 + (l + /3 2 )(l + g 2 



6V2ip 



2^ + (3 



■\l+D(#,P)]. 
(126) 

It can be readily verified that in the Hall limit these re- 
sults reduce to the corresponding expressions presented in 
Section 4.3. 



5.4 Dissipation rate 

The results of Section 4.4, given by equations (121) 
and (122), continue to apply also in the Ohm regime. Note 



11 This condition implies a minimum inclination angle of the sur- 
face poloidal field to the rotation axis of ~ 49°, as compared with 
30° for the ideal MHD-like case (b rh > 1/V3). 



12 We have implicitly assumed that the Hall parameters arc > 
in this case; this is explicitly justified in Section 6. 
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that the second of the above equations can be written in the 
form 



(3 ■ e ')o = -^2 (l + a/3 2 ) [1 + D(^,(3)} . (127) 



5.5 Results 

The results of the foregoing analysis are collected in Ta- 
bles 3 and 4, which are patterned on Tables 1 and 2, re- 
spectively, in Section 4 (where the corresponding data for 
the Hall regime are presented) . The implications of these re- 
sults are discussed in conjunction with those of our findings 
for the Hall and ambipolar-diffusion regimes in Section 6. 



5.6 Formulation in terms of scalar conductivity 

As we noted in Section 2.1, the analysis of the problem in 
the Ohm regime can be alternatively carried out in terms 
of the scalar conductivity ao (equation 11). In the case of a 
two-component plasma with q <C 1 this is just the "electron" 
electrical conductivity a c (equation 27). With our adopted 
normalization (see equation 66), the corresponding dimen- 
sionless conductivity at the midplane of the disc is 



O~c0 



47t/itc s ct c hTC s 



dAo 



a 



(128) 



Ao 



From equation (128) it is seen that, in the Ohm regime, 
= a ( 2 )Cr . 13 This makes it possible to relate the results 
derived in this section to the equivalent scalar-conductivity 
formulation. In particular, Cases (i) and (ii) identified above 
correspond to a\a a being > 1/2 and < 1/2, respectively. 
Resistive-MHD disc models that have appeared in the liter- 
ature typically use Ohm's law in the form E — —v X B/c + 
J/a e and ignore the second (Hall) term on the right-hand 
side of the more general expression given by equation (25). 
On the face of it, this might be justified by the fact that the 
ratio of the conductivity prefactors of the second and first 
terms on the right-hand side of this equation, | <th | co/ax , 
is ~ \Pc\, which is <JC 1 in the Ohm regime. However, this 
argument does not account for the fact that the Hall term in 
the expression for E' is perpendicular to J, which implies, 
in particular, that the Hall term appears in the prefactor of 
the radial current density J r in the normalized expression 
for the azimuthal electric field component E$ (equation 94) . 
The only other contributor to this prefactor at z = is as- 
sociated with the advective (ccdxB /c) term, which could 
in principle be subordinate to the Hall-current term. Now, 
these two terms (Hall and advective) in equation (94) corre- 
spond to the two terms (/? and 2tp, respectively) in the nu- 
merator of the expression for \db r /db^\o (equation 115) that 
are used in our solution classification scheme for the Ohm 
regime. In Cases (i) and (ii) Yo|Ao| = ao<?eo/|/?eo| > 1/2 
and so the Hall term can be neglected, but in Case (Hi) this 
inequality is reversed and the Hall term dominates. 



This also follows directly from the relation (87) in this limit. 



6 DISCUSSION 

The first term within the first pair of parentheses on the 
right-hand side of equation (122) represents the ohmic con- 
tribution to the Joule dissipation in the disc, whereas the 
second term corresponds to the ambipolar dissipation. (The 
Hall term in Ohm's law, i.e., the term oc J X B in equa- 
tions 25 and 26, does not contribute to the J • E' dissi- 
pation.) This suggests (see Konigl 1997) that the slip fac- 
tor so = l/q/3 2 — AoAo can be used to distinguish be- 
tween ambipolar diffusion-like (AoAo > 1) an d Ohm-like 
(AoAo < 1) sub- regimes in the nominal Hall parameter do- 
main (|Ao| < 1 < |/3 c o|). 14 For comparison, both AoAo > 1 
and |Ao| > 1 must hold in the actual ambipolar diffusion 
regime, whereas both AoAo < 1 and \/3 e o\ < 1 must be sat- 
isfied in the genuine Ohm regime. This criterion forms the 
basis of the classification given in the second column of Ta- 
ble 1. The relevance of this criterion is demonstrated by the 
fact that the first row of inequalities (Case i) exactly repro- 
duces the corresponding constraints in the ambipolar diffu- 
sion case (see WK93 and Konigl 1997) and that the third 
row (Case Hi) is identical to the first row of inequalities in 
Table 3 (Case i in the Ohm regime). 15 

The third column in Table 1 classifies the solutions 
according to the value of the magnetic coupling parame- 
ter Ao (the Elsasser number, which in the Hall regime is 
equal to Yo|Ao|)- Although the constraints obtained for 
large (Ao > 1/2; Cases i and Hi) and small (Ao < 1/2; 
Cases ii and iv) values of Ao are clearly distinct in their 
forms, it is interesting to note that the sub-Keplerian ro- 
tation and wind-launching requirements (the first and sec- 
ond inequalities, respectively) together imply that To > 1 
in all cases. This inequality already follows directly from 
the conditions T ]Ao| > 1/2 and |Ao| < 1 that define the 
Hall Cases (i) and (Hi), but it is noteworthy that all viable 
solutions in the Hall domain satisfy this constraint, which 
was previously inferred in the ambipolar diffusion regime 
(see WK93). This result can be visualized by plotting the 
sub-Keplerian-rotation and wind-launching constraints for 
each of the Hall sub-regimes. First, however, we present the 
relevant curves in the ambipolar diffusion limit (Fig. 1; cf. 
Fig. 3 in WK93). The dark, solid lines in this figure show 
the run of b r versus bj,, computed from equation (106), for 
To = 10, 2 and 1. In fully sub-Keplerian discs db r /db<p must 
be < below the base of the wind (i.e. for z < z\,\ see 
equations 105 and 106). This condition constrains viable so- 
lutions to lie to the left of the long-dashed curve that marks 
the locus of points where dbr/db,/, = 0. On the other hand, 
these solutions must also satisfy the wind- launching crite- 
rion (equation 111), which constrains them to reside above 
the light solid line for the given value of To. It is evident 
from these plots that no physically viable solutions exist for 
To < 1. The corresponding curves for Cases (i) — (iv) in 
the Hall limit are shown in the four panels of Fig. 2: they 
confirm that viable solutions satisfy To > 1 also in all the 
Hall sub-regimes. 

As was shown in Section 5, the magnitude of the pa- 



14 This can also be seen directly from equation (34), which im- 
plies that s = rjA I no ■ 

15 These correspondences apply also to the solution characteris- 
tics listed in Tables 2 and 4. 
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Table 3. Parameter constraints for wind-driving disc solutions in the limit where the Ohm diffusivity dominates. Three distinct cases 
can be identified, as indicated. The meaning of the four inequalities in each sub-regime is the same as in Table 1. The parameters tp, (3 
and qf3 used in the text are related to the parameters employed in the table through tp = To, f3 = 1/fto an d iP = l//?e0- 
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Table 4. Key properties of viable disc solutions in the Ohm regime. See the caption to Table 2 for a description of the listed quantities. 
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rameter combination T |/3h>| compared to 1/2 is one of the 
classification criteria also in the Ohm regime (see the second 
column in Table 3). When To|/3io| > 1/2 the second classifi- 
cation criterion in the Ohm regime can again be expressed in 
terms of the Elsasser number Ao (which in this limit equals 
To/3eo/3io) being either > 1/2 (the Ohm Case i) or < 1/2 
(the Ohm Case ii; see the third column in Table 3). Simi- 
larly to the situation in the Hall Cases (i) and (Hi), the fact 
that |/3io| < 1 implies that To must be > 1 also in these two 
Ohm sub-regimes. The third Ohm sub-regime is defined by 
To|Ao| < 1/2 and T jA=o| > 1. Since |/3 c0 | < 1 in the Ohm 
domain, it follows that To > 1 in this case too. 

The parameter To (which is equal to Ao in the ambipo- 
lar diffusion regime) thus turns out to be of fundamental 
importance in the theory of diffusive wind-driving discs in 
that it is inferred to be > 1 for all viable solutions irre- 
spective of the conductivity regime. Physically, the condi- 
tion To > 1 expresses the requirement that the momentum 
exchange time of the neutrals with the particles that dom- 
inate the momentum of the charged species (the compar- 
atively massive "ions" in our formulation) be shorter than 
the dynamical time (i.e. the Keplerian orbital time). This re- 
quirement is evidently more basic than having the nominal 
neutral-field coupling parameter Ao be > 1. In fact, as noted 
above, some of the distinct parameter sub-regimes that we 
identified are defined by the inequality A < 1/2. However, 
even in the latter cases Ao remains bounded from below by 
a numerical factor that typically is not <JC 1 (Ao > «o/3 
and Ao > ao/V& in the Hall and Ohm regimes, respectively, 
from the wind launching condition; see Tables 1 and 3). 

The lower bounds on Ao reflect our assumption that the 
modelled discs are strongly coupled, which implies, in par- 
ticular, that the transverse magnetic field component starts 
to grow already at z — (see the solution curves presented 
in Paper II). As was first demonstrated by Li (1996) and 
Wardle (1997), viable solutions of discs in which the bulk 
of the matter is not well coupled to the field can also be 



constructed. In these solutions, significant field-line bend- 
ing that enables a centrifugally driven wind to be launched 
from the disc surface only starts well above the midplane. 
Although the basic disc properties above the height where 
b r starts to increase rapidly are similar to those of strongly 
coupled configurations, this location depends on the de- 
tailed density and ionization structure of the disc and can- 
not be readily determined a priori. It roughly coincides with 
the height where \db r /db^\ (estimated using equation 115), 
which is typically <C 1 at the midplane under these circum- 
stances, has grown to a value ~ 1. However, disc material 
already moves inward below this height on account of the 
fact that l&^l becomes < 1 when b r is still <C 1, which means 
that a significant magnetic torque (of density oc B^dB^/dz) 
can be exerted even at lower elevations. 

The magnitude of the field-strength parameter a is 
bounded from below by the sub-Keplerian rotation con- 
straint (the first inequality in Tables 1 and 3) and from above 
by the wind-launching constraint (the second inequality in 
Tables 1 and 3). If the magnetic field is too weak it will 
not drive an outflow but instead mediate a "two-channel" 
flow within the disc (see Salmeron et al. 2007). On the other 
hand, if the field is too strong it will not bend out sufficiently 
to drive a centrifugal wind. Our typical solutions are charac- 
terized by a < 1. We have found that the ratio |& r b/&^>b| of 
the radial to azimuthal magnetic field amplitudes at the disc 
surface is > 1 in the ambipolar diffusion regime (WK93) as 
well as in all the Hall sub-regimes (Section 4.3). However, 
in the Ohm regime this is true only for Case (i), whereas in 
the other two Ohm sub-regimes \b$b\ > b r b (see Table 4). 
Nevertheless, even in the latter two cases b r b is required to 
exceed a certain lower bound (which is at least 2/\/3; see 
Section 5.2) to satisfy the wind launching condition. It is 



16 Such configurations were later found to arise naturally in cer- 
tain models of disc formation from the collapse of rotating molec- 
ular cloud cores (see Krasnopolsky & Konigl 2002). 
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Figure 1. Plots of b r versus b$ (each normalized by (2T0) 1 / 2 ) 
in the quasi-hydrostatic region of the disc (dark solid lines) for 
To = 1, 2 and 10 in the limit where ambipolar diffusion dom- 
inates. Physically viable solutions must satisfy two constraints. 
First, they must lie to the left of the long-dashed line, which 
marks the point where dbr/db^ changes sign: this ensures that 
the azimuthal velocity is sub-Kcplcrian within the disc (i.e. w^, 
given by equation 105, is < for z < Second, they must also 
lie above the light solid line for the corresponding value of To , so 
that the wind-launching criterion (equation 111) is satisfied. It is 
evident from the figure that these two constraints together imply 
To > 1. 



also worth noting in this connection that, even if the mag- 
netic field at the disc surface is not bent strongly enough to 
satisfy the wind launching requirement, vertical magnetic 
angular-momentum transport along a large-scale field that 
threads the disc could potentially still take place by other 
means, such as magnetic braking (see footnote 1). 

By combining the wind- launching and wind-loading 
constraints (the second and third inequalities, respectively, 
in Tables 1 and 3) and using the expression for h in the 
fifth column of Tables 2 and 4, one finds that h < 1, i.e. 
that the vertical confinement of the disc is primarily mag- 
netic (due to the vertical gradient of (B 2 + B^)/8tt) rather 
than tidal for all the cases in the Hall and Ohm regimes. 
This result, which was previously shown by WK93 to apply 
to ambipolar diffusion-dominated discs, is corroborated by 
the full solutions presented in Paper II. Predominantly mag- 
netic confinement is thus seen to be a generic property of our 
wind-driving disc model. The generality of this result was, 
however, questioned by Ferreira (1997), who suggested that 
it might be a consequence of the simplification pv z =const 
adopted in our treatment of the mass conservation relation. 
In addressing this question it is interesting to note that the 
inequality h < 1 was inferred in this paper already on the 
basis of the hydrostatic approximation, which does not make 
use of the continuity equation. On the other hand, our ar- 
gument involved an estimate of z b , the normalized height of 



the disc surface, which corresponds to a disc region where 
the hydrostatic approximation may no longer be accurate. 
A definitive resolution of this issue would require using a 
more elaborate model in which the above simplification is 
not made. 

The wind-loading and energy dissipation constraints 
(the third and fourth inequalities, respectively, in Tables 1 
and 3) provide lower and upper bounds, respectively, on the 
radial velocity parameter e. The lower bound reflects the fact 
that the height of the disc surface scales as h whereas the 
scale height scales as 1/e, so that z b /h, which is required 
to be > 1, scales as e 2 (see the fifth and sixth columns in 
Tables 2 and 4). The upper bound is a consequence of the 
fact that both the electric field and the current density scale 
as e (see the seventh column in Tables 2 and 4), so that the 
Joule dissipation (which scales as e 2 ) might exceed the rate 
of gravitational energy release (oc e) if the inflow speed were 
too high. The solutions presented in Paper II have values of 
e in the range ~ 0.3 — 1. With these parameters one infers 
a nominal emptying time < 10 4 yr for a representative pro- 
tostellar disc, which is consistent with the mean duration of 
the earliest (Class 0) phase of protostellar evolution but is 
very much shorter than the total age of accreting protostars. 
This is the essence of the "short evolution time" critique of 
pure-wind angular-momentum transport models of proto- 
stellar discs, which has led to the argument that alternative 
angular-momentum transport mechanisms must dominate 
in such systems (e.g. Shu et al. 2008). 17 Note, however, that 
even if a centrifugally driven wind does not control the over- 
all evolution of the disc throughout its lifetime, it might still 
dominate the angular momentum transport in regions of the 
disc from which strong outflows are launched, which may be 
especially pertinent during the FU Orionis-type rapid ac- 
cretion episodes that are thought to account for most of the 
mass deposition on to low-mass protostars (e.g. Calvet et 
al. 2000). Furthermore, wind- mediated angular momentum 
transport could in principle also contribute in regions where 
a radial transport mechanism dominates (see Salmeron et al. 
2007). It is also worth noting that significantly lower mean 
inflow speeds could be attained in more realistic, vertically 
stratified wind-driving disc models (e.g. Konigl & Salmeron 
2009), and in particular in weakly coupled discs (see Li 1996 
and Wardle 1997). Further work is, however, required to de- 
termine the extent to which this effect would increase the 
predicted disc evolution time. 

We have been careful to use the absolute values of $0 
and /3e0 in delineating the various sub-regimes in the Hall 
and Ohm domains in view of the fact that our definition 
(equation 4) of the Hall parameters has an explicit depen- 
dence on the magnetic field polarity and thus allows them to 
have both positive and negative values. However, in our ta- 
ble entries for the Hall Cases (ii) and (iv) and for the Ohm 
Case (Hi) we implicitly treated the ion and electron Hall 
parameters as positive quantities. This can be readily jus- 
tified by observing that these cases correspond to \/3\ being 
> 2ip and taking note of the inequality /? > — 2ip (equa- 



17 These alternative mechanisms potentially include vertical 
transport along a large-scale open magnetic field by means other 
than the steady wind that has been the focus of our discussion; 
see footnote 1. 
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Figure 2. Same as Fig. 1, but in the limit where the Hall diffusivity dominates. The four panels correspond to Cases (i) through (id) 
in Table 1. Note that, as in the ambipolar diffusion regime shown in Fig. 1, physically viable solutions satisfy To > 1 also in each of the 
Hall sub-regimes. 



tion 108) that all viable solutions must obey. These two 
inequalities can be simultaneously satisfied only if /3 > 0. 
In other words, for these parameter sub-regimes a self- 
consistent wind-driving disc solution can be obtained only 
if the field has a positive polarity but not if sgn{B z } < 0. 
The sensitivity to the field polarity reflects the explicit de- 
pendence on the Hall conductivity in these cases (see Sec- 
tions 2.1 and 5.6) and is a specific prediction of our model 
that is checked in the numerical work presented in Paper II. 
We emphasize that our discussion here concerns the nature 
of certain sub-regimes, in which only solutions with (3 > 
are predicted to exist. This property is distinct from the 
dependence of wind-driving wind solutions on the field po- 
larity in other parameter sub-regimes in which both positive 
and negative values of /3 are allowed but in which chang- 
ing the sign of f3 modifies the behaviour of the solution (see 
WK93). An example of the latter dependence is given in Pa- 
per II for the Hall Case {in), which admits f3 values in the 
range —tp/2 < (5 < 2ip (see equation 114). 

Our analysis has concentrated on discs with a charged- 



particle composition that can be represented as a two- 
component plasma consisting of electrons and comparatively 
massive ions, with the parameter q (equation 7) being <g; 1. 
This composition is appropriate to relatively low-density 
disc regions or even to denser zones that do not contain 
dust (see footonote 5). Weakly ionized disc regions that have 
densities > 10~ 13 gcm~ 3 and contain dust may be better 
approximated by a two- component, q = 1 plasma corre- 
sponding to oppositely charged grains of equal mass (see 
Section 2.1). In this case the Hall regime is not present (see 
equation 95) and the parameter space of viable solutions 
consists of only two regimes (ambipolar diffusion and Ohm) . 
One can readily verify that, under these circumstances, the 
Ohm domain contains only two sub-regimes, the analogues 
of the Ohm Cases (i) and (ii) for a q <JC 1 plasma. 18 This 



18 Note from equations (95)- (97) that the parameter ip = To 
always appears in the conductivity tensor in the combination (1 + 
q) tp. Therefore, the parameter Tq in the expressions for the g«l 
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finding is not surprising in view of the fact (see Section 5.6) 
that the Ohm Case (iii) for a g < 1 plasma arises from 
the Hall-current term in Ohm's law, which is missing when 
9 = 1. 

The viability constraints that we derived on diffusive ac- 
cretion discs in which the entire excess angular momentum is 
transported by a centrifugally driven wind can be useful also 
for generalizing the model to incorporate other modes of an- 
gular momentum transport. In particular, one can argue (see 
Salmeron et al. 2007) that radial angular-momentum trans- 
port through MRI-induced turbulence could operate at disc 
locations where the sub-Keplerian rotation constraint (ex- 
pressed by the first inequality in Tables 1 and 3, in which 
the various parameters are taken to have their local, density- 
dependent values) is violated. For example, in the ambipolar 
diffusion case (as well as in the Hall Cases i and ii and the 
Ohm Cases i and ii) this corresponds to 2Ta 2 oc p\/ p being 
< 1. As this quantity is expected to increase with height 
on account of the increase in the fractional ionization of the 
gas on going away from the disc midplane, one can envi- 
sion a situation in which 2Ta 2 is < 1 (with both radial and 
vertical angular-momentum transport taking place) in the 
vicinity of the midplane but becomes > 1 (with only ver- 
tical transfer into a wind remaining relevant) closer to the 
disc surfaces. This criterion can be employed for construct- 
ing a "hybrid" disc model in which both vertical and radial 
angular momentum transport take place (through a large- 
scale, ordered magnetic field and a small-scale, disordered 
field, respectively; see Salmeron et al. 2007 for details). 



7 CONCLUSION 

We have investigated a protostellar accretion disc model in 
which the dominant angular momentum transport mecha- 
nism is a centrifugally driven wind launched along a large- 
scale, ordered magnetic field that threads the disc. The ba- 
sic assumptions we adopted are the same as those of the 
radially localized disc model constructed by WK93, and our 
main purpose has been to extend the analysis presented in 
that paper, which focussed on discs in which ambipolar dif- 
fusion dominates the magnetic diffusivity inside the disc, 
to a general diffusivity regime. In particular, we included 
also the Hall and Ohm diffusivities, which are relevant to 
higher-density disc regions than those covered by the WK93 
treatment. We employed a tensor-conductivity scheme that, 
in conjunction with an ionization-balance calculation, can 
be used to determine realistic vertical conductivity profiles 
for protostellar discs in the context of this model. However, 
in this paper and its follow-up, Paper II, we adopt a simpler 
treatment and assume that the three model parameters that 
are needed to identify the distinct solution regimes for an 
ion-electron plasma remain constant with height in the disc. 
The three variables employed for this purpose in the current 
paper are the electron and ion Hall parameters (/3 e and 
respectively, with \f3i\ <C |/? e |) and the neutral-ion coupling 
parameter T (equation 82). In Paper II we use instead two 
independent ratios of the conductivity tensor components 

Cases (i) and (ii) in Tables 3 and 4 is replaced by 2To in the 
corresponding expressions for q = 1. 



as well as the Elsasser number A (see items (Hi) and (iv) in 
Section 3.13). 

We determined the parameter regimes where physically 
viable disc solutions could be found by employing the hy- 
drostatic approximation and imposing the following require- 
ments (originally applied in the ambipolar diffusion regime 
by WK93 and Konigl 1997) on "strongly coupled" systems 
(discs in which the magnetic field is dynamically well cou- 
pled to the bulk of the matter already at the midplane): (i) 
the disc (in contradistinction to the wind) rotates at sub- 
Keplerian speeds; (ii) a wind launching condition (which 
yields a lower limit on B r /B z ) is satisfied at the disc surface; 
(iii) most of the disc material is located below the wind- 
launching region; and (iv) the rate of Joule dissipation does 
not exceed the rate of gravitational energy release in the 
disc. Our main results can be summarized as follows: 

• In addition to the ambipolar diffusion regime considered 
in WK93, there are four distinct sub-regimes in the Hall 
diffusion-dominated parameter domain and three distinct 
sub-regimes in the Ohm domain. (A fourth potential Ohm 
sub-regime was found not to be consistent with the wind- 
launching requirement.) The four Hall sub-regimes naturally 
divide into two ambipolar diffusion-like and two Ohm-like 
parameter regions: In the former pair one of the sub-regimes 
has the same structural properties as the ambipolar diffusion 
regime, whereas in the latter pair the solution characteris- 
tics of one of the sub-regimes are identical to those of one 
of the Ohm sub-regimes. In the case (not treated in detail 
in this paper) of a high-density disc plasma dominated by 
oppositely charged grains of equal mass, only an ambipolar 
diffusion regime and an Ohm regime (with two sub-regimes 
corresponding to our Ohm Cases i and ii) are present. 

• All viable solutions, irrespective of the diffusivity 
regime, are found to satisfy T > 1. The physical requirement 
expressed by this condition, that the neutral-ion momentum 
exchange time be shorter than the disc orbital time, is thus 
indicated to be a fundamental constraint on wind-driving 
discs of this type. 

• Viable solutions are also characterized by the parame- 
ters ao (midplane Alfven-to-sound speed ratio) and e (mid- 
plane inflow-to-sound speed ratio) being < 1. However, 
weakly coupled discs, characterized by values of A (the 
midplane value of the Elsasser number; see Section 3.13) 
that are <C 1, could have ao <C 1 and e <C 1. In such sys- 
tems a significant bending of the magnetic field away from 
the rotation axis (required for driving a centrifugally driven 
wind) takes place only above some finite height in the disc, 
although vertical angular momentum transport along the 
large-scale field typically occurs already at lower elevations. 
Although our analysis does not directly apply to weakly cou- 
pled discs, our results should be useful for understanding the 
basic qualitative aspects of their behaviour. 

• The wind-launching and wind-loading conditions (the 
aforementioned requirements (ii) and (iii), respectively) to- 
gether imply that magnetic squeezing (by the gradient of 
the magnetic pressure force associated with the transverse 
field components B r and B$) dominates the gravitational 
tidal compression of the disc. A more elaborate model is, 
however, needed to verify the full generality of this result. 

• The transverse magnetic field at the disc surface is 
dominated by the radial component in all the Hall param- 
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eter sub-regimes (including their "twins" in the ambipolar- 
diffusion and Ohm domains), and by the azimuthal com- 
ponent in the remaining two Ohm sub-regimes. However, 
even in the latter two cases the poloidal surface magnetic 
field must be sufficiently strongly bent to satisfy the wind 
launching requirement. 



Centrifugally driven winds are not the only means of 
vertical angular momentum transport in a protostellar disc 
(magnetic braking is an example of another option), and 
there are also alternative mechanisms that involve radial 
transport along the plane of the disc. It is nevertheless 
instructive to investigate the limiting case in which wind 
transport dominates at a given radius in view of the ubiq- 
uity of energetic outflows in protostellar systems and the 
inferred association of disc winds with FU Orionis-type 
high- accretion-rate events. The constraints considered in 
our model can also be useful for identifying the parame- 
ter regimes where other mechanisms (such as MRI-induced 
turbulence) are likely to operate, possibly even at the same 
radial location as wind-mediated angular momentum trans- 
port (Salmeron et al. 2007). 

The parameter regimes that we identified were deter- 
mined under the assumption that the values of the param- 
eters used in the classification scheme did not vary with 
height in the disc. This approximation becomes less accurate 
as one moves from the ambipolar-diffusion regime to the Hall 
regime and then (at even higher disc densities and column 
densities) to the Ohm regime. In particular, in the latter 
case it can be expected that if a real disc could be described 
as being in this regime at the midplane then it would transit 
to the Hall regime further up and would likely be ambipolar 
diffusion-dominated at the surface. For this reason (and also 
because wind-driving discs might not be massive enough to 
harbour an Ohm regime within their weakly ionized zones) 
we have emphasized the Hall regime in our discussion and 
do not consider the Ohm domain in Paper II. However, our 
findings for the Ohm regime could still be applicable to the 
innermost, collisionally ionized regions of protostellar discs 
(and possibly also to disc regions participating in an FU 
Orionis-type outburst), in which an anomalous ohmic resis- 
tivity might dominate. 

The actual parameter regime that characterizes a given 
section of an accretion disc of the type that we have mod- 
elled will be determined by the global magnetic flux distri- 
bution, the density structure (which depends in part on the 
mass accretion rate) and the ionization profile (which de- 
pends in part on the nature of the ionization sources as well 
as on the disc composition). It is conceivable that some of 
the sub-regimes that we have identified are not commonly 
realized in Nature or that they correspond to unstable con- 
figurations. More insight into these questions would likely 
emerge from comparisons with observations as well as from 
additional analytic and numerical work. Our study has, how- 
ever, indicated that, if such wind-driving accretion flows are 
indeed present in real protostellar systems, they will exhibit 
certain generic properties (including good neutral-ion cou- 
pling, subthermal Alfven and inflow speeds and magnetic 
squeezing within the field-line bending region) irrespective of 
the particular parameter regime to which they correspond. 
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APPENDIX A: THE e B « APPROXIMATION 

In their consideration of a weakly ionized disc containing 
a two-component plasma, WK93 noted the qualitative sim- 
ilarity of solutions characterized by the same value of the 
parameter combination (e — cb). They attributed this result 
to the fact that the only change in the underlying system of 
equations introduced by switching to a reference frame that 
moves with the radial velocity VBrO = — £bc s of the poloidal 
flux surfaces involves the magnitude of the torque that is 
required to remove the excess angular momentum, so that 
all radial-velocity terms (except in the angular momentum 
equation, which remains unchanged) are modified simply by 
subtracting VBrO from the radial velocity component. We 
appealed to this result in Section 3.13 to motivate setting 
£b = in our analysis. 

One could however question this approach on the basis 
of the following consideration. Starting from Ohm's law in 
the form E' = E + v X B/c = (Att/c 2 )t] ■ J, where rj is 
the diffusivity tensor, we rewrite it in terms of normalized 
quantities as 



WBrO — w r 



m -fjp \ I 0.5 W r 

fjp ??H / V 2U>0 



, (Al) 



where fjp = fjo + VA and where 770, Vh and f)A are ob- 
tained by multiplying the expressions in equations (29)-(31), 
respectively, by 1/tiTCsa 2 ,. We have also written w Br o = 

fBro/cs = WE4.0 and Wb^O = r(Q,B0 — 51k)/c s = — WErO 

for the normalized midplane radial and azimuthal veloci- 
ties, respectively, of the poloidal flux surfaces. By inspect- 
ing equation (Al) one may be induced to infer that, since 
the off-diagonal terms of the diffusivity tensor dominate in 
the ambipolar-diffusion and Ohm regimes whereas the diag- 
onal terms are dominant in the Hall limit, one should retain 
the latter terms when considering the Hall regime and the 
former terms when treating the other limits. This inference 
could be carried over to the analysis of equation (98) in the 
text, in which the matrix is related to the inverse of the ma- 
trix that appears in equation (Al). Specifically, specializing 
to the midplane (where b r — — 0) and assuming q <C 1 for 
simplicity, one can use equations (95)-(97) to deduce that 
Ai — —A4 oc ap/a± (with the factor qfi 2 in equations 99 
and 102 representing the ratio 770 /??a). One can similarly 
deduce that the term f3 in the expressions for A2 and A3 
is oc <th/oj_ (with the term involving rp in equations 100 
and 101 arising from the advective contribution to the elec- 
tric field). Adopting the above reasoning, it would appear 
that one should retain the diagonal terms A\ and A4 in the 
ambipolar-diffusion and Ohm regimes and the off-diagonal 
terms A2 and A3 in the Hall limit. However, this procedure 
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cannot be implemented if one sets ee to be identically zero in 
equation (98), as we have done in our analysis, which casts 
doubt on the validity of the approach followed in the text. 

In addressing this question it is important to keep in 
mind the basic attributes of the system that we wish to 
model. We are primarily interested in the quasi-steady be- 
haviour of a disc in which mass and poloidal magnetic flux 
that were originally part of a large-scale equilibrium con- 
figuration of a molecular cloud core are carried in by the 
accretion flow. In this picture, the accretion is enabled in 
large measure by the vertical magnetic transport of angu- 
lar momentum, and the inward motion of the matter is, in 
turn, responsible for the radial transport of the field lines 
(which, however, generally lag behind the matter because 
of the disc's diffusivity) . The vertical magnetic transport of 
angular momentum implies that J r B z > and hence that 
w r < (see equation 92), whereas the inward bending of 
the magnetic field lines by the inflowing gas implies that 
J<f>B z > and hence that < (see equation 91). The 
fact that the motion of the field lines is controlled by that 
of the gas implies that kjb^> = -t»Er would also be < 
(see Section 3.12). In a similar vein, we expect the radial 
velocity of the poloidal field lines (described by the param- 
eter WBrO = — £b) to satisfy w r o < WBrO (or, equivalcntly, 
£B < £). 19 

We now proceed to examine the consequences of the 
alternative formulation outlined above in light of the just- 
discussed constraints. To this end, we substitute for j r o and 
j^o in equation (98) from equations (92) and (91), respec- 
tively, to obtain 

Al TOErO + A 2 e B = ^ , (A2) 
Zip 

A 3 w Er « + A 4 e B = -^^ . (A3) 

V 

In the ambipolar-diffusion and Ohm regimes we retain only 
the A\ and A4 terms. In the ambipolar diffusion limit 
Ai = — A4 ~ 1 and K « 1 + ip 2 ~ ip 2 (since ip is equal 
in this case to the midplane value of the neutral-field cou- 
pling parameter A, which is expected to be > 1), so equa- 
tion (A3) implies eb = < 0, which is inconsistent 
with our expectation that e_g > 0. Similarly, in the Ohm 
limit Ai = — A4 « q/3 2 and K « q/3 2 + ip 2 , so equation (A3) 
implies eb = 2(ip/q[3 2 + qp) 2 /ip)w</,o, which is again inferred 
to be < 0, contrary to our physical expectation. In a similar 
vein, we retain only the A2 and A3 terms in the Hall regime. 
Using equations (100) and (101) and focussing on the case 
q <§C 1, we infer from equations (A2) and (A3) that 




and 



WEK) 




One can readily verify that equations (A4) and (A5) can- 
not be simultaneously satisfied while also respecting the 
constraints < eb/e < 1 and wet-o/w^o < 0. These ar- 
guments indicate that a formulation based simply on the 
dominant terms in the conductivity tensor is not physically 
self-consistent for the problem under consideration. 

In contrast, one can verify that the formulation adopted 
in the text is self-consistent - in the sense that the neglected 
terms in equations (A2) and (A3) never come to dominate 
the ones that were retained - under the assumption eb /e < 
1 . In particular, using the expression (113) for WErO and 
equations (99)-(102), it is straightforward to ascertain that 
the ratios |A2eB/^4iWEro| and |AieB/^4.3WEro| remain < 1 
in this case for all the viable-solution parameter regimes 
identified in Sections 4 and 5. 

This paper has been typeset from a TfjjX/ I^TffjX file prepared 
by the author. 



19 It is conceivable that €3 and e could have opposite signs, but 
such a situation would likely only arise in the context of a lo- 
calized, time-dependent phenomenon. An example of this pos- 
sibility is provided by the cyclical behaviour exhibited by the 
core-collapse solutions of Tassis & Mouschovias (2005) near the 
boundary of their computational "central sink." This behaviour 
is indeed localized (in both space and time) and clearly does not 
correspond to the global, quasi-steady accretion process analyzed 
in this paper. 



